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PURIFICATION OF SUBSTANCES BY SLOW FRACTIONAL 
FREEZING 


Frank W. Schwab and Edward Wichers 





ABSTRACT 


This paper describes two technics for purifying substances by slow fractional 
freezing. 

The first involves the slow lowering of a cylindrical cell, filled with the fused 

substance, through a heating coil in such a way that freezing begins at the bottom 
and progresses upward, as the cell emerges from the coil, until the whole mass has 
solidified. The solidified column, tube and all, is cut into the desired fractions. 
The second technic, which is suitable for larger quantities, affords better 
control of the rate of freezing and provides a larger solid-liquid surface for a given 
quantity of material. A spherical flask filled with the liquid is cooled at a con- 
trolled rate in such a way that freezing begins at the wall of the flask and proceeds 
inward at a regulated slow rate until the desired fraction has solidified, after 
which the remaining liquid is siphoned out of the flask. 
The purification of benzoic acid was much more rapid by this method than in 
earlier trials by crystallization from solvents. Acetanilide was also effectively 
purified. Observations incidental to the work show that the freezing point of 
pure acetanilide is 114.29 +0.01° C. 



























In an earlier paper ' brief reference was made to the use of fractional 
freezing as a means of purifying benzoic acid. This work was done 
oa very small scale. It followed closely a technic previously used 
by R. S. Jessup, of the Heat Measurements Section of this Bureau. 
A cylindrical tube, 2 cm in diameter, containing 50 g (45 to 50 ml) 
of the fused acid (mp 122.4° C) was slowly lowered through a 
heating coil. As the tube emerged from the coil, freezing began at 
the bottom and progressed until the whole mass was frozen. During 
the freezing, which extended over a period of 6 hours, the liquid 
portion was constantly stirred. 

The foregoing procedure was also used to purify benzoic acid on a 
larger scale, for quantities of about 500 g. The glass tubes used were 
46 cm in diameter, and the column of molten acid was about 30 cm 
long. The tube was lowered through the heating coil by supporting 
iton a platform floating in a tank of water, as shown in figure 1. 
Water was allowed to flow slowly from the tank through an outlet 
valve on a lower floor of the building. This was done to keep the 
hydraulic head substantially constant. Some difficulty was experi- 
enced with stoppage of the needle valve by particles of dirt. Alto- 





'F. W. Schwab and E. Wichers, Preparation of benzoic acid of high purity, J. Research NBS 25, 747 (1940 
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gether, there were enough minor problems in operating this mechani- 
cally simple device to suggest that a somewhat more elaborate arrange. 
ment, such as a clock- or motor-driven movement, would have proved 
more satisfactory. 

The liquid portion of the acid was stirred with a stream of nitrogen 
gas flowing from a Branham-Sperling bubbler.? The total period of 
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Ficure 1.—Sketch of apparatus for the purification of a substance by slow fractional 
freezing in a cylindrical cell 


A, Heating element backed by a layer of insulating material supported on an 7 Pyrex cylinder; B, 8 
Branham-Sperling bubbler tip for stirring with nitrogen; C, cylindrical cell; D, cork support for cell; 
E, float guide; F, float; G, supports; H, container for water; and I, water outlet. 


freezing was from 16 to 20 hours. An advantage of this method is 
that the acid is readily separated into fractions by allowing all of it to 
freeze and then cutting the column, tube and ail, into suitable lengths. 

A quantity of benzoic acid, which had an initial purity of 99.91 
mole percent, was fractionally frozen three times in this apparatus, 
and each time the last one-fourth frozen was discarded. The yield of 
purified acid thus was about 40 percent. Its purity was 99.997 mole 


2J, R. Branham and E. O. Sperling, Bubbler tip of Pyrex glass for difficult absorptions, J. Research NBS %, 
701 (1989) RP1214. v . - su 
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percent, as determined by measuring its freezing range in the manner 
described in an earlier publication.® 

In the study of a process of this sort, it is helpful to know what 
fraction of the total impurity is removed by a single operation, in this 
case, @ single crystallization. This fraction will be referred to as the 
efficiency of the process of purification. The over-all increase in 
purity can be used to calculate the efficiency of the process of purifica- 
tion. In this instance the substance contained, in 100,000 moles of 
starting material, 90 moles of impurities susceptible of removal by 
fractional freezing. After the third step of fractionation, the content 
of impurities was 3 moles in 100,000, or 3.3 percent of the amount 
originally present. Assuming equal efficiency of separation in all 
three steps of purification, the fraction of impurity not removed in a 
single step is the cube root of 0.033; that is, 0.82. The average effi- 
ciency of purification in this instance was therefore 68 percent. 

When benzoic acid was fractionally frozen in the manner described, 
the solid material was clear and did not show evidence of entrapment 
of liquid or other macroscopic contamination. It seems probable 
that the relatively low efficiency of purification was caused by failure 
to maintain equilibrium at the solid-liquid interface. This could re- 
sult from irregularities in the seemingly smooth progress of the cell 
through the heating coil, from irregular cooling caused by drafts or 
temperature changes around the apparatus, or simply from too rapid 
freezing. As freezing takes place the impurities increase in concen- 
tration in the very thin layer of liquid immediately adjacent to the 
solid surface. It is well known that this layer is not mixed with the 
body of the liquid by stirring, hence impurities can be transferred 
from it to the body of the liquid only by diffusion. It seems evident 
therefore that the degree of approach to ideal separation of impurities 
must be a function of the rate of deposition of crystalline material per 
unit area of surface. 

Instead of attempting to improve the relative efficiency of purifica- 
tion by freezing in cylindrical cells, another apparatus was prepared 
in which the area of freezing surface was larger in relation to the 
quantity of material to be frozen, and the rate of freezing could be 
more steadily maintained and more easily regulated. In this appa- 
ratus it was possible also to handle much larger quantities than can 
well be done in cylindrical vessels. The apparatus is shown in figures 
2and 3. It consisted essentially of a long-necked, 5-liter flask of 
Pyrex glass imbedded in an especially effective thermal insulating 
material, silica aerogel, containing about 15 percent of silicon,‘ and 
provided with a small resistance heater wound around the base of 
the neck. The flask was filled with molten benzoic acid at a tempera- 
ture about 10° C above its freezing point, enough to allow the loss of 
heat to the surroundings to become steady before freezing began. 
The liquid was agitated with a glass stirrer. The stem of the stirrer 
passed through ‘a hole in the flat bottom of a glass cylinder which 
closely fitted the neck of the flask, and which was placed with its 
'¥. W. Schwab and Edward Wichers, Precise Measurement of the Freezing Range as a Means of De- 
termining the Purity of a Substance, Temperature—Its Measurement and Control in Science and Industry, 
p. 256-264 (Reinhold Publishing Corporation, New York, N. Y., 1941). 

The impurities measured by the freezing range are limited to those soluble in the liquid but insoluble in 
the crystals. The same limitation applies, of course, to a process of purification based on fractional freezing. 
‘Ind. Eng. Chem. 81, 827 (1939). 

' The glass stirrer is made by sealing two rectangular flats about 1 cm by 6cm toa glassrod. The long 
edges of the flats are centered and sealed to opposite sides of the rod, forming an angle at 20 to 25 degrees 


with it. A stirrer of this type vibrates less than the usual propeller type. The authors are indebted to 
W. G. Schlecht for the design of this stirrer. 
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bottom near the upper end of the resistance heater. The bottom of 
the sleeve formed a baffle to prevent loss of heat by convection in the 
neck of the flask and also sublimation of the acid into the neck. Ax 
the flask was imbedded in insulating material to a point above the 
heater on the neck, heat could be added to the system at this point 
to regulate to a limited extent the net loss of heat from the body of 
the flask. With the heater furnishing 5.8 watts, the rate of total loss 
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Fieurs 2.—Apparatus for the purification of a substance by slow fractional freezing 
in a round-botiom flask. 


A, glass stirrer; B, sleeve to prevent loss of heat by convection; O, five-liter flask; D, heater for balancing 
the loss of heat through the neck of the flask; E, insulating material (silica aerogel); F, container for in- 
sulating material; and G, substance being purified. 


of heat from the flask was about 7.5 watts, as determined from the 
amount of benzoic acid frozen and the latent heat of fusion. The 
rate of loss of heat was also roughly determined by the rate at which 
the liquid benzoic acid cooled just before it began to freeze. 

After the desired fraction of the acid was frozen, the stirring as- 
sembly was removed and the remaining liquid was drawn off by suc- 
tion. Then, a point-source radiant heater in the form of a small 
45-watt lamp (automobile headlight type), with its electrical leads 
enclosed in a glass tube, was lowered through the neck to the center 
of the flask. Since benzoic acid is nearly opaque to red and infrared 
radiation, the heat from the lamp caused the surface layer of the acid 
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to melt. About 15 minutes’ exposure to the lamp was required to 
melt 100 g. The surface of the solid which had been in contact with 
the impure mother liquor was thus washed by drawing off the remelted 
portion. This washing was done three times. The solid acid formed 
a nearly uniform layer about 2.8 cm thick over the whole inside of the 
flask except that it was thinner near the neck, where the level of the 
liquid had subsided because of the reduction in volume which occurred 
during freezing. The acid was nearly transparent except for shrink- 
age cracks which developed on cooling to room temperature. 











FiacurE 3.—Assembly for partial melting and washing of crystallized material. 
H, Radiant heater; I, layer of crystals. 


In the foregoing manner about 10 kg of benzoic acid was twice 
frozen, each time rejecting about 40 percent of the quantity in the 
flask. ‘The combined products of the first fractionation of two 5-kg 
lots thus yielded more than enough to fill the flask for the second 
fractionation. The average purity of the product of the first frac- 
tionation was 99.986 mole percent and that of the second 99.998 mole 
percent. Since the starting material had a purity of 99.91 mole per- 
cent, it will be seen that the relative efficiency of separation was about 
85 percent in each step. This figure is comparable with the efficiency 
of 68 percent reported for the cylindrical cell, because nearly the same 
fraction was rejected. It represents a marked gain in efficiency. 
This is attributable in part, perhaps, to a steadier rate of freezing, but 
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mainly to a slower rate of freezing. In the flask this was approxi- 
mately 0.18 g/cm? of crystallizing surface per hour (a rate of growth 
of 0.024 mm per minute), whereas in the cylindrical cell it was 1.25 
g/cm? of crystallizing surface per hour. It is obvious that, for any 
given quantity of liquid to be frozen, the use of the whole surface of a 
spherical container as a solid-liquid interface offers a marked advan- 
tage over freezing at one end surface in a cylindrical vessel. It is true 
that as freezing progresses in the flask the surface exposed to the liquid 
becomes smaller, thereby increasing the rate of freezing per unit area, 
However, the rate of loss of heat can be decreased to balance this change, 
or even to over balance it, so as to offset the increasing concentration of 
impurities in the liquid. Crystalline substances in general are such 
good conductors of heat relative to materials used for thermal insulation 
that the growing layer of solid does not greatly decrease the rate of loss 
of heat at the wall of the flask. All of these effects are decreased as 
the size of the flask is increased. Within the range of available sizes of 
flasks, the principal limitation to the scale on which this method of 
purification can be used is the problem of getting the frozen material 
out of the flask. This can be done by heating the whole flask in an air 
or liquid bath, or by the use of a solvent. An internal heater similar 
to that described, but of larger capacity, could also be used. In the 
present instance the thick layer of benzoic acid was melted by care- 
fully warming the flask with a gas flame. If this is done, care must 
be taken to keep the material that first melts from freezing again in 
the space between the solid and glass, otherwise subsequent thermal 
expansion of the solid cake will break the flask. 

A comparison between the efficiencies of purification here obtained 
by fractional freezing with those obtained in earlier work (footnote 1) 
on crystallization from solvents may be of interest. In that work the 
average efficiency of purification in a single crystallization from water 
was 14 percent and from benzene, 29 percent. ‘These efficiencies 
were computed on the basis that in a series of n steps of purification 
the mean fraction of impurity not removed in a single step is the 
nth root of the fraction of the original impurity that remains after 
steps. With these efficiencies 11 recrystallizations from benzene, in 
the manner there followed, or 25 from water, would have been needed 
to achieve the same increase in purity that was accomplished by the 
two fractional freezings. In considering this comparison it should be 
borne in mind that the purity of the starting material used in the 
recrystallization from solvents was considerably higher than that 
used for the fractional freezing and that the character of the impurities 
was not necessarily the same. It should also be noted that the effi- 
ciencies in question do not necessarily represent the highest that could 
be obtained by the respective methods. 

The apparatus used for benzoic acid was also used (by W. H. 
Smith and W. W. Walton, of this Bureau) to purify a quantity of 
acetanilide. A single fractionation, in which about 50 percent was 
rejected, resulted in an increase in purity from 99.80 mole percent 
(determined from its freezing range) to 99.98 mole percent. The 
relative efficiency of purification therefore was 90 percent. In this 
instance the solid acetanilide was removed from the flask by dissolving 
it in hot benzene. On cooling the benzene solution the acetanilide 
was recovered in the desired form of small crystals. The indicated 
purity of the product therefore involved one step of crystallization 
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fom benzene in addition to the fractional freezing. It was de- 
termined after removing the benzene trapped by the crystals. The 
extent of this incidental contamination by the solvent amounted to 
0.11 percent for the crystals as dried at room temperature and 0.05 
percent after the crystals were crushed somewhat. The amount of 
solvent in the substance was determined by a method, developed by 
the authors ® for determining water and other volatile material in 
pure substances. The freezing temperature of this product, after 
removing benzene from it, was 114.28° C. The purity of the sub- 
stance and the accuracy of the measurement of the temperature were 
suficient to warrant giving 114.29°+0.01° C as the freezing tempera- 
ture of pure acetanilide. 

It is believed that slow fractional freezing deserves more considera- 
tion as a means of purification than it’ has received heretofore for sub- 
stances stable at temperatures somewhat above their freezing points. 
The slow and steady rate of crystal growth that can be maintained 
oflers a marked advantage over crystallization from solvents, as 
commonly practiced, with respect to approaching equilibrium between 
the solid and liquid phases. The method need not be limited to 
substances which freeze above room temperature. It should be 
equally applicable, with a suitable refrigerant, to substances which 
freeze at lower temperatures. 


The authors express their appreciation to C. L. Gordon for the 
preparation of the drawings appearing in this paper. 


Wasnineron, April 1, 1944. 


‘4 sample of acetanilide was fused in a glass apparatus, which was evacuated through a trap immersed in 
liquid air, The trap was then closed off from the rest of the apparatus, and the liquid air was removed, 
allowing the trap to warm up to room temperature. The amount of benzene removed from the acetanilide 
was estimated from the pressure, measured by means of a mercury manometer, and the temperature and 
volume of the closed system. 
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ASTUDY OF THE PROPERTIES OF HOUSEHOLD BLANKETS! 


By Herbert F. Schiefer, Hazel Tharp Stevens,? Pauline Beery Mack,® and 
‘ Paul M. Boyland 





ABSTRACT 


The fiber composition, weight, thickness, compressibility, compressional resil- 
jence, thermal transmission, air permeability, breaking strength, and shrinkage of 
156 different blankets are recorded. The effects of laundering; of laundering and 
renapping; of laundering, renapping, and abrasion; of dry cleaning and renappng; 
and of dry cleaning, renapping, and abrasion on these properties of a large number 
of blankets are shown. A linear relationship was found between the compres- 
sional resilience and the wool content of cotton-wool blankets. The thermal 
transmission of the blankets was found to be independent of the kind of fiber. 
The reciprocal of thermal transmission was found to be related linearly to the 
thickness. The thermal transmission computed by means of the equation 
1/T=3,0 t&.; +0.63, where 7' is the thermal transmission in Btu/(°F hr ft?), 
and t,, is the thickness in inches at a pressure of 0.10 lb/in.?, was found to agree 
with the measured values within +10 percent, 95 times out of 100. Empirical 
relationships were also found among thermal transmission, thickness at 1.0 
lb/in.4, and we. ressibility ; and among -thickness at 0.10 lb/in.*, compressibility, 
~ weight. e relation between breaking strength and weight, and that be- 
tween Nonkiag strength and compressibility, are discussed. Minimum require- 
ments are suggested for the properties of blankets for use in a performance 
specification. 
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I. INTRODUCTION 


Blankets are an essential article in every home, hospital, and hotel, 
ard in other places of abode. In time of war they form a vital part 
of the protective clothing of the soldier, sailor, and marine on battle 
fronts and in training camps. They furnish shelter to the wounded 
and to those forced to take to lifeboats and life rafts at sea. In time 
of peace blankets are an essential article for camping and other out- 
door recreation and sport. In a major emergency they provide 
temporary shelter to the unfortunate victims. 

Although blankets form such an essential part of life in war and in 
peace, adequate methods for testing them and essential information 
regarding their performance characteristics have not been available. 
New blankets are usually soft, lofty, fluffy, flexible, and warm. It 
is desirable to have them retain these characteristics during use. 
Changes in the characteristics of blankets during use are probably 
attributable mainly to the effect of laundering or dry cleaning. This 
is particularly true if the temperature and alkalinity of the water 
used in laundering wool or acetate blankets are too high or not con- 
trolled, or if the blankets receive excessive mechanical action during 
washing or dry cleaning. 

New blankets differ considerably in their characteristics, and some 
will change more than others when washed or dry cleaned. A blanket 
which is considerably felted when new, for example, would not be 
expected to change greatly, whereas one which is soft, fluffy, and flexi- 
ble is more likely to become felted and less flexible during cleaning. 
The cleaning procedure for blankets of the latter type would have to 
be controlled more closely to minimize changes in their desirable 
characteristics. 

The results of tests on a large number of blankets of different fiber 
composition are given in this paper. The effects of laundering; of 
laundering and napping; of laundering, napping, and abrasion; of 
dry cleaning and napping; and of dry cleaning, napping, and abrasion 
on the properties of blankets are discussed. Information is given 
concerning the correlation between various properties of the blankets, 
and requirements for a performance specification of blankets are 
suggested. 

II. MATERIALS TESTED 


The blankets which were tested, 156 in number, are listed in table 
1 according to their fiber composition and weight. Most of the 
blankets were obtained directly from manufacturers by the National 
Bureau of Standards, the Medical Corps of the United States Army, 
and the Ellen H. Richards Institute, Pennsylvania State College. 
Some were purchased from retail stores. A few were submitted for 
test by Government departments. It is believed that these blankets 
are a good representation of the qualities of blankets available at the 
time this investigation was made. 











Results of testis on 1456 blankets 


TABLE i.— 
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III. TESTING PROCEDURE 


The weight, thickness, compressibility, compressional resilience, 
breaking strength, shrinkage, air permeability, and thermal trans- 
mission of the blankets before and after washing or dry cleaning were 
measured with equipment* and methods of testing® previously 
described. 

Weight is expressed in ounces per square yard hereafter referred 
to as areal weight. Thickness is given for stated pressures. Breaking 
strength is expressed as the load at failure in the “grab” procedure, 
wherein a specimen considerably wider than the clamps is tested. 
Shrinkage is expressed as the percentage change in the dimensions 
of the specimen in each of its two principal directions. 

Compressibility, (At/t,.9)/AP, is the decrease in thickness, Af, in unit 
thickness at a pressure of 1.0 lb/in.’, t,o, per unit change of pressure, 
AP. For convenience in the numerical calculation of compressibility, 
Atis taken to bet) ;—t, 5, where f).,and tf, 5.are the thickness values at pres- 
sures of 0.5 and 1.5 lb/in.’, respectively, and AP, the difference between 
these two pressures, is therefore 1.0 lb/in?. The unabridged unit, 
(in./in.)/(Ib/in.?), rather than the contracted form, in. ?/Ib, previously 
used, is used in this paper because its physical meaning is more 
readily understood. 

Compressional resilience is the work recovered when the pressure 
is decreased from 2.0 to 0.10 lb/in? expressed as a percentage of the 
work done when the pressure is increased from 0.10 to 2.0 lb/in?. 

Thermal transmission is measured by the rate of loss of heat energy 
through a unit area of the material per unit temperature difference, 
when there is a temperature difference of 20°C (68°F) between the hot 
plate and the hood (cold plate). 

Air permeability is measured by the volume of air flowing in unit 
time through a unit area of the material when the pressure difference 
across one thickness of the blanket is 4 in. of water. 

The blankets were conditioned by exposure to air having a relative 
humidity of 65 percent and a temperature of 70°F, with a tolerance 
of +2 percent in relative humidity and +2°F in temperature. All 
= tests, except thermal transmission, were made under these con- 

itions. 

Some of the blankets were washed 10 times in a commercial laundry 
without brushing or renapping after washing. Some were washed 
10 times in the United States Army laundry at Fort Myer, Va., 
without renapping after washing. Others were washed 10 times at 
the National Bureau of Standards in a reversing wash wheel, in an 
0.5-percent soap solution at 100° F for 15 min, followed by three 
5-min rinses also at 100° F. The samples were then centrifuged for 
1 min and dried in a horizontal position without tension. Some of 
these blankets were not renapped after each washing, whereas others 
were renapped after each washing. Duplicate samples of the latter 
were renapped and then abraded Sees the next washing to simulate 
the mechanical wear during use. Some of the blankets were dry 


* Herbert F. Schiefer, The Compressometer, an instrument for evalucting the thickness, See, and 
compressional resilsence of textiles and similar materials, BS J. Research 10, 705 (1933) RP561. 
P A zs ae Riots improved apparatus for measuring the thermal transmission of teztiles, J. Research NBS 
Herbert F. Schiefer and Paul M. Boyland, Improved instrument for measuring the air permeability of fabrics, 
J. Research NBS 28, 637 (1942) RP1471. 
' Herbert F. Schiefer 


, Advantages of a blanket-and-sheet combination for outdoor use, J. Research NBS 30, 
208 (1943) RP1529. 
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ceaned in a commercial plant with a washer 36 in. in diameter and 
4 in. in length and normal load of 60 to 70 lb. Stoddard solvent, 
conforming to Commercial Standard Specification CS3-40, and a 
commercial paste type of soap, which contained about 10 percent of 
water, were used. One pound of soap was used for every 25 lb of 
lad. The dry cleaning procedure consisted of ® 10-min run in 
slvent without soap; a 20-min run in solvent and soap; and a 20-min 
rinse in solvent. The blankets were centrifuged and then dried in a 
tumbler for 30 min. The temperature of the air leaving the tumbler 
was approximately 160° F. The dry-cleaned samples of blankets 
were renapped after each dry cleaning. Duplicate samples were 
rnapped and then abraded before the next dry cleaning. 

The properties of the blankets were measured when new and most 
of them were measured again after 1, 5, and 10 washings or dry 
cleanings. 

Some of the laundered and dry-cleaned blanket samples were 
abraded with a machine similar to the one developed by Haven.° 
The Haven machine was modified to permit testing a blanket sample 
14 in. wide and 18 in. long. The length of the stroke was reduced 
014 in. The area abraded was 14 by 14 in. The blanket sample 
was fastened to two clamps specially constructed for these tests. 
One clamp was fastened to the reciprocating mechanism of the 
machine. The blanket sample was placed over a horizontal, smooth 
shaft, which was 1 in. in diameter and mounted in ball bearings, and 
to the other clamp was attached a 5-lb weight. A horizontal bar 
weighing 5 lb 7 and covered with No. 8 duck was placed on the blanket 
directly above the shaft. This bar was rectangular in cross section 
and was held in place by vertical guides at each end. The blanket 
sample was drawn back and forth over the shaft, which rotated 
freely, and underneath the duck-covered bar, which remained sta- 
tionary. The back-and-forth movement of the blanket sample 
relative to the duck fabric abraded the nap of the blanket. After 
200 back-and-forth movements, the blanket sample was turned over 
and the reverse side was abraded the same number of strokes. This 
abrasion procedure was repeated before each subsequent laundering 
or dry cleaning. 

The fiber composition of the blankets was determined microscopi- 
cally. Ifmore than one kind of fiber was found in a blanket, a chemical 
analysis was made. 


IV. RESULTS AND DISCUSSION 


The results of the tests of the original blankets are given in table 1. 
The results of the tests of the blankets after laundering and dry clean- 
ing are omitted to conserve space; however, the over-all changes in 
the properties of the blankets are brought out in the following 


discussion. 
1. EFFECT OF LAUNDERING 


The effect of laundering without renapping on the properties of 43 
blankets is shown in figure 1, where the amount of change in a property 
after 1, 5, and 10 washings, expressed as a percentage of the value of 
the property of the blanket when new, is plotted as abscissa. For 


‘Textile World 76, p. 2654-2656, 2662 (1929). 
’In some of the early tests a weight of 1 Ib was used but was found to be insufficient. 
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FIGURE 1.—Frequency with which the blankets changed a given amount in properties after 1, 5, and 10 washings. 
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each property the cumulative frequency, expressed as the percentage 
of blankets that changed more or less than a given amount is plotted 
as ordinate, the left-hand frequency scale being used for positive 
changes and the right-hand scale for negative changes. As an ex- 
ample, in 10 washings the areal weight decreased slightly for less than 
\0 percent of the blankets, whereas half of them increased 15 percent 
or more in areal weight, and a sixth of them increased in areal weight 
as much as 25 percent. 

All of the blankets shrank in the warp and filling directions when 
washed. Most of the shrinkage occurred in the first washing, although 
the shrinkage increased consistently with the number of washings. 
The shrinkage in the warp direction was approximately twice as great 
as the shrinkage in the filling direction. The shrinkage of the blankets 
accounts for their increase in areal weight. It also accounts in part 
for the increase in breaking strength and thickness and for the decrease 
inair permeability. In addition to shrinkage, the blankets became 
felted or matted during washing. This condition is measured by the 
change in compressibility. A high value for compressibility indicates 
a greater amount of napping, whereas a low value indicates a greater 
amount of felting or matting. The considerable decrease in com- 
pressibility with the number of washings indicates that the blankets 
became felted or matted. This change affects the feel to the hands 
and also the appearance of the blanket. Brushing or renapping of 
the blanket after washing decreased these changes and restored the 
blanket more nearly to its condition when new. A blanket which 
became felted during washing increased in breaking strength and 
decreased in air permeability. 

The compressional resilience of blankets after the fifth washing is 
appreciably lower than after the first or tenth washing. A high value 
for compressional resilience indicates greater ability of the blanket to 
come back to its initial state upon release of a compressive load. The 
compressional resilience of a blanket will be decreased if the fibers be- 
come damaged during use or laundering. This decrease results be- 
cause damaged fibers have less ability to recover from flexure than 
undamaged fibers. On the other hand, the compressional resilience 
of a blanket will increase if it becomes felted during use or laundering. 
This increase results because there is less motion of the fibers relative 
toeach other when a compressive load is applied owing to the greatly 
reduced compressibility of a felted blanket. The energy which is 
lost, frictional forces between fibers times relative movement of fibers, 
is therefore considerably reduced and the compressional resilience, 
which is the ratio of the energy recovered when the compressive load 
is removed to the total energy expended when the compressive load 
is applied, is increased. The consistent increase in compressional 
resilience between the fifth and tenth washings resulted from felting 
ofthe blankets. The discussion of the change in resilience of a blanket 
during laundering is general and independent of the kind of fiber. 

he compressional resilience of a new blanket depends, however, 
upon the kind of fiber from which it is made. This dependence of 
resilience upon the fiber composition will be discussed later in this 
paper. 

_ When a blanket shrinks the thickness at a pressure of 0.10 Ib/in.? is 
Iereased, and when it becomes felted the density is increased and the 
thickness at 0.10 lb/in.* is decreased. However, if further shrinkage 
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occurs during felting the thickness may be increased. As all of they 
changes occurred when the blankets were washed, the thickness 4 
0.10 lb/In.? showed no consistent increase or decrease. This is no} 
true for the thickness at a pressure of 1.0 lb/in.*, because the thicknes 
at this pressure is not greatly affected by felting. Laundering ip. 
creased the thickness at 1.0 lb/in.*, most of the increase being obtained 
in the first washing owing to the shrinkage. 

The thermal transmission of the blankets was affected little by the 
number of launderings. None of the blankets increased over 10 per. 
cent in thermal transmission, and less than 10 percent of the blankets 
decreased over 10 percent during the 10 washings. A low value fo 
thermal transmission indicates a high insulating value. In gener! 
the thermal transmission of a blanket changes inversely with both 
thickness and compressibility. Hence if thickness increases whil 
compressibility decreases, the changes tend to annul each other in 
their effect on thermal transmission. This probably accounts fo 
the fact that the thermal transmission of the blankets is in generil 
affected little by laundering. 

The large increase in the filling breaking strength for some 10 per. 
cent of the blankets was obtained because they had a low strength 
when new. These blankets were highly napped and had a high 
compressibility. It is noteworthy that sos initially weak blankets 
did not become weaker when laundered but became stronger. This 
increase resulted from the felting and shrinkage produced by launder. 
ing. 

It is apparent that a number of properties of the blankets changed 
materially during laundering. Some of the changes, such as the 


increase in breaking strength, are not objectionable. The smallness 
of the change in thermal transmission is particularly desirable. Hov- 
ever, the rather large decrease in compressibility is objectionable. 
It denotes a change in the feel and appearance of the blankets. Fortv- 
nately, most of this change can be overcome by brushing or renapping 
of the blanket after washing. This fact is brought out in subsection 2. 


2. EFFECT OF LAUNDERING AND RENAPPING, AND OF LAUNDER. 
ING, RENAPPING, AND ABRASION 


The effects of laundering and renapping; and of laundering, renap- 
ping, and abrasion on the properties of 27 blankets are shown in 
figure 2, where the amount of change in a property after 10 launder- 
ings, expressed as a percentage of the value of the property of the 
blanket when new, is plotted as abscissa, and the cumulative fre- 
quency is plotted as ordinate as in figure 1. The changes in a real 
weight, thickness, thermal transmission, air permeability, breaking 
strength, and shrinkage are similar in general to the changes showi 
in figure 1 for laundering without renapping. Compressibility 1 
materially improved by renapping after laundering, whereas compres- 
sional resilience is decreased by renapping. The effect of renapping 
after washing on the appearance of several typical blankets is shows 
in figures 3 and 4. 

The curves showing the changes in the properties of the laundered, 
renapped, and abraded blankets differ only slightly from the curves 
for the unabraded blankets. The difference in appearance between the 
abraded and unabraded blankets was also slight. The weight and 
thickness of the abraded blanket samples are consistently lower, and 
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Ficure 2.—Frequency with which the blankets changed a given amount in properties afler 10 washings and 10 nappings; 
10 washings, 10 nappings, and 9 abrasions;10 dry cleanings and 10 nappings; and 10 dry cleanings, 10 nappings, 


and 9 abrasions. 
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FIGURE 3.—Appearance of blankets when new and after the tenth washing before and 
after the tenth napping. 


Left column, original. Top row, 100 percent of wool. 
Middle colum, 10 washings and 9 nappings. Middle row, 51 percent of wool, 49 percent of cotton. 
Right column, 10 washings and 10 nappings. Bottom row, 100 percent of cotton. 
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Figure 4.—Appearance of blankets when new and after the tenth washirg before and 
after the tenth napping. 


Left column, original. Top row, 23 percent of wool, 31 percent of cotton, 49 
percent of viscose rayon. 

Middle column, 10 washings and 9 nappings. Middle row, 6 percent of wool, 62 percent of cotton, 32 
percent of viscose rayon. 

Right column, 10 washings and 10 nappings. Bottom row, 12 percent of wool, i percent of cotton, 
87 percent of viscose rayon. 
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the thermal transmission and air permeability are consistently higher 
than for the unabraded samples. 


3. EFFECT OF DRY CLEANING AND RENAPPING, AND OF DRY 
CLEANING, RENAPPING, AND ABRASION 


The effects of dry cleaning and renapping, and of dry cleaning, 
renapping, and abrasion on the properties of 27 blankets are also 
shown by curves in figure 2. These curves show practically the same 
changes as those for the laundered blankets. 


4, EFFECT OF FIBER .COMPOSITION ON THE PROPERTIES OF 
BLANKETS 


The blankets in table 1 are listed according to fiber composition. It 
can be readily seen that there is no correlation between the properties 
of the blankets and the fiber composition, except for compressional 
resilience. The average compressional resilience, computed for those 
blankets having similar fiber composition, is plotted in figure 5 against 
the wool content. A linear relationship is obtained for blankets con- 
taining a mixture of wool and cotton. The difference in compressional 
resilience between the groups of blankets is statistically highly signi- 
ficant, except between the 5-percent-wool and the all-cotton groups. 
The addition of 5 percent of wool does not increase the compressional 
resilience significantly. If a portion of the cotton is replaced with 
viscose rayon, the compressional resilience is significantly decreased, 
and if a portion of the cotton is replaced with acetate rayon, the 
compressional resilience is increased slightly. Several of the 100- 
percent-wool blankets were known to be made from all-virgin wool. 
The average compressional resilience of this group was 50 percent, 
which is higher than the general average of the all-wool blankets. 
Likewise, a number of blankets were known to be made from a high 
percentage of reprocessed and reused wool. The average compres- 
sional resilience of this group was only 41 percent, which is considera- 
bly below the general average of the all-wool blankets. There seems 
to be a relationship between the compressional resilience of blankets 
and the quality of wool used. If mechanically damaged wool is used 
in blankets, the compressional resilience of the blanket is lowered. 


5. RELATIONSHIP BETWEEN THERMAL TRANSMISSION AND 
THICKNESS 


The reciprocal of the thermal transmission of blankets and some 
underwear fabrics recently measured in this laboratory is plotted in 
figure 6 against the thickness at a pressure of 0.10 Ib/in?. A linear 
relationship is obtained. Lines indicating 5 and 10 percent deviation 
are shown in figure 6. The relationship between thermal transmis- 
sion and thickness is given by equation 1 

, 1 

7 =3.0%,,+-0.03 bul CF hr ft*), (1) 
where T is the thermal transmission, and t&, is the thickness at a 
pressure of 0.10 lb/in?. It will be noted that 7 is 1.6 when &,; is 
zero. This is the value which is obtained when a measurement is 
made with the bare apparatus. This value of 7 was used in com- 
puting the thermal insulating values given in table 1. The thermal 
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Figure 5.—Linear relationship between compressional resilience and wool content of 
blankets made from miztures of wool and cotton. 


Each point represents the average compressional resilience of blankets having approximately the same 
fiber mixture. 
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transmission computed by means of equation 1 is within +10 percent 
of the measured value about 95 times out of 100. This difference is 
frequently less than the difference which is obtained when two dif- 
ferent portions of one blanket are measured. 

The relationship between the thermal transmission and the thick- 
ness at a pressure of 1.0 lb/in. ? is given in equation 2. 


1 


T= 75t,,40.71 


Btu/(F hr fe). (2) 


This equation yields values for thermal transmission which are within 
+10 percent of the measured value about 75 times out of 100. The 
residuals between observed and computed values are related to the 
compressibility. The equation can therefore be modified to give 
equation 3, 

1 


T= 7544-071 


—0.650+0.20, (3) 


including a term for the compressibility, C, of the material. Equa- 
tion 3 also yields a value of 1.6 for 7’ when the values for t,o and C 
are zero, and corresponds to the measurement made with the bare 
apparatus. The values of 7’ computed from equation 3 are within 
+10 percent of the measured values about 92 times out of 100. 

Equations 1, 2, and 3 are applicable only when the thickness is 
measured at the specified pressures. However, formulas could be 
derived for thicknesses measured at other pressures. The equations 
are applicable only when the thermal transmission is measured with 
the equipment previously described. If the thermal transmission 
values are measured with some other equipment, then it is possible 
to compare these measured values with those that would be obtained 
with the equipment described in Research Paper RP1055, as computed 
by either equation 1 or 3. 

Equation 1 may be written as equation 4, 


R=3.0t)1-+0.63(°F hr ft?)/Btu, (4) 


where R is the reciprocal of 7 and may be defined as thermal resistance. 
The constant 0.63 is a thermal resistance which is equivalent to an 
effective thickness of 0.21 in. and is caused by the layer of still air at 
the exposed surface of the specimen. This value is comparable with 
the value of about 0.16 in. reported by Baxter and Cassie.® Equation 
4 shows that two layers of a textile material do not yield twice the in- 
sulation of one layer of the material. Tests were made on a number of 
blankets by using one and then two layers. The results agreed very 
well with the allies computed from equation 4 or 1. 

It is apparent from the rather large constant in equation 4, equiva- 
lent to an effective thickness of 0.21 in. for to.,, that it is impractical to 
attempt to measure the effect of construction and other factors on the 
thermal resistance of thin materials, such as light underwear and 
dress fabrics, since the thermal resistance contributed by thin fabrics 
may amount to only some 10 to 20 percent of that of the layer of 
still air. 


*R. 8. Cleveland, J. Research NBS 19, 675 (1937) RP 1055. 
® mal insulating properties of clothing, J. Textile Inst. 84, No. 7, T41-T'54 (July 1943). 
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6. RELATIONSHIP BETWEEN WEIGHT AND THICKNESS 


A linear relationship between thickness and areal weight is often 
assumed or taken for granted. It is shown in figure 7, however, that, 
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FicurgE 7.—Relationship between the thickness of blankets at 0.10-lb/in.? pressure 
and the weight. 


such a correlation is very poor. The best straight line through the 
plotted points corresponds to equation 5 


ty 1=0.0155W, (5) 


in which ¢, ; is the thickness at 0.10 lb/in.?, and W is the areal weight 
A further study indicates that the points are scattered about this 
line in a manner closely associated with the compressibility. If the 
ratio of the observed values of t)., to the values computed by means of 
equation 5 are plotted against compressibility, an approximate linear 
relationship is obtained. Equation 5 can thus be modified to yield 
equation 6, 


ty 1. =0.05CW, (6) 


where C is the compressibility. 
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If the thicknesses of the blankets computed by means of equation 6 
are plotted against the measured values, figure 8, a straight line with 
unit slope is obtained. The points, although still scattered consider- 
ably, lie for the most part within +20 percent of the line given by 
equation 6. An approximate value of the weight, W, can thus be 
computed if the values of f, and C are known. Equation 6 is useful 
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The dashed lines indicate a 20-percent departure. 








as a guide in specifying minimum compatible requirements for a 
performance specification for blankets. 


7. RELATION OF STRENGTH TO WEIGHT AND TO COMPRESSIBILITY 


The breaking strength of blankets was found to increase in general 
with the 8 weight and to decrease with the compressibility. 
Although the functional relationships between strength and areal 
weight and between strength and compressibility are not well defined, 
they suffice for use as a guide in specifying minimum requirements for 
a performance specification of blankets that will be compatible with 
the other requirements and also provide reasonable assurance of 
satisfactory performance. 
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8. RELATIVE THERMAL INSULATION PER UNIT WEIGHT 


The reciprocal of thermal transmission has been used as a measure 
of the relative thermal insulation of textiles. This value divided by the 
areal weight of the material has been used to indicate the efficacy as 
an insulating material and is given in table 1. It can be readily seen 
that it varies with the thickness but not with fiber composition. In 
general, the greater the thickness of the blanket the lower is the value 
of the thermal resistance per unit weight. The correlation, however, 
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Figure 9.—Linear relationship between the reciprocal of thermal transmission times 
weight and the value computed by means of equation 7. 


is far from good for the blankets as a group. The correlation is 
greatly improved if the blankets are grouped according to their com- 
pressibility. The general relation can be expressed by equation 7, 
where 7’ is thermal transmission, W is weight, C is compressibility, 


p= 0-40—0Aty 1 +.0.05, (7) 


and t), is thickness at 0.10 lb/in?. The values of 1/7W computed by 
means of equation 7 are plotted in figure 9 against the values computed 
from the measured values of 7’ and W. The points lie fairly close to a 
straight line having a slope of unity. This indicates that equation 7 
represents the data fairly well. 
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The physical significance represented by equation 7 is not imme. 
diately apparent, except that the relative thermal insulation of 
blankets per unit weight depends upon the thickness and compres. 
sibility and therefore upon the structure of the blanket. However, 
the compressibility and thickness of a blanket of a given weight ar 
inversely related to the density. The relative thermal insulation 
per unit weight given by equation 7 is therefore inversely related to 
the density. This relationship can be derived directly from equa- 
tion 1 or 4 and is given by equation 8. 


1 _ 3%. ; 0.63 
TW Ww 8 


The ratio W/t); is equal to 126, where 6 is the density of a blanket in 
pounds per cubic foot when it is under a pressure of 0.10 Ib/in®. By 
substituting 126 for W/t)1, equation 9 results. 


1 0.25 , 0.63 


TW 3° WwW’ 0) 


A good approximation of 1/7W can be computed from equations 7, 
8, and 9. Equation 7 requires only measurements of thickness and 
compressibility, which can be made on the finished blanket in a few 
minutes and without cutting and destroying the blanket. The value 
computed by means of equation 8 or 9 is a better approximation, but 
requires determining the areal weight in addition to the thickness at 
0.10 Ib/in?. These equations have a practical application in the 
establishment of minimum compatible requirements for the proper- 
ties of blankets in a performance specification. 


9. THERMAL INSULATING VALUE 


The thermal insulating values, J, given in table 1 were computed by 
means of equation 10 


I=100,— (10) 


where 7’ is the thermal transmission of a blanket, and 7’, is the value 
obtained with the bare apparatus, namely, 1.60 Btu/(°F hr ft*). By 
substituting 1.60 for 7, and the value of 7 from equation 1 in equation 
10, equation 11 is obtained, 


I 100 ae 


~ 1+0.21/to.1’ 

where ¢, is the thickness at a pressure of 0.10 lb/in?. In figure 10 
the values of J from table 1 are plotted against t, and the curve 
shown is the locus of equation 11. The thermal insulating value 
obtained for underwear fabrics, blankets, and several lined jackets 
agree very well with the curve representing equation 11. It can be 
readily seen that the thermal insulating value of textiles having the 
same thickness, t) ;, as measured in this investigation, does not depend 
upon the fiber composition and type of fabric construction. A very 
good value of J can be computed with equation 11 if the value of 
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t; is known. This value of a fabric or garment is readily measur- 
able and does not require cutting or destroying the material for the 
test. 
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FiaurEe 10.—Relationship between the thermal insulating value, I, of blankets and 
the thickness, ty, at 0.10-lb/in® pressure. 


The curve represents equation 11. 


V. SUGGESTED REQUIREMENTS FOR A PERFORMANCE 
SPECIFICATION FOR BLANKETS 


The results obtained in this investigation can be used as the basis 
for the preparation of performance specifications for blankets. Limits 
for compressional resilience, weight, compressibility, thickness at 
0.10 lb/in.? pressure, breaking strength, and shrinkage are suggested 
in the following paragraphs. 


1. COMPRESSIONAL RESILIENCE 


A soft, fluffy, flexible, and warm blanket having a high compres- 
sional resilience is expected to retain these characteristics more nearly 
during use than one having a low compressional resilience. The com- 
pressional resilience was found to be linearly related to the fiber compo- 
sition of blankets containing cotton and wool. It could therefore 
serve as a requirement in a performance specification in place of fiber 
composition, since it is readily measured without cutting or destroying 
a blanket. The minimum values suggested for a performance speci- 
fication are given in table 2. The approximate fiber composition 
corresponding to the suggested values for compressional resilience is 
indicated. The suggested values for compressional resilience are 
applicable to blankets regardless of their weight. 
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TABLE 2.—Suggested values of compressional resilience for use in a performance 
specification 








Approximate fiber 
Minimum composition 
com pres- 
Type of blanket sional nies l 
resilience ool, 
minimum | Cotton 








Percent Percent Percent 

45 90 10 
40 35 
35 60 
30 85 
25 100 




















. WEIGHT 


The weight of a blanket fabric depends primarily upon the conditions 
prevailing for the use of the blanket. For purposes of specification, it 
is suggested that the minimum values of the weight of the fabric in 
ounces per square yard be the even numbers, starting with 4. 


3. COMPRESSIBILITY 


The compressibility of a blanket is a measure of its loftiness or degree 
of napping. The values to be specified depend again upon the condi- 
tions of intended use and also upon individual preferences for the feel 
and appearance of blankets. For purposes of specification, it is sug- 
gested that the minimum values specified be in increments of 0.05 
beginning with 0.20 (in./in.)/([b/in?). A blanket having a low com- 

ressibility is one which is greatly felted, little napped, and stiff, or 
sone ay and one having a high compressibility is lofty, soft, highly 
napped, and flexible. 
4. THICKNESS 


The thickness of a blanket is probably the most important character- 
istic to be specified in a performance specification. It has been shown 
that the thermal transmission and the thermal insulationg value cor- 
relate exceedingly well with the thickness at 0.10 lb/in?. Simple 
formulas, equations 1 and 11, have been derived for calculating these 
values from this thickness. Minimum values for thickness are sug- 
gested in table 3 for various weights and compressibilities in preference 
to either thermal transmission or thermal insulating values, because 
of the ease of measurement of thickness which is obtained when the 
measurements for compressional resilience and compressibility are 
made, and because the thickness can be measured without the need 
of cutting or destroying the blanket. However, computed thermal 
transmission and thermal insulating values are given in table 4 for 
different thicknesses as a guide for converting the thickness values 
of table 3 to these quantities. The densities in pounds per cubic foot 
at a pressure of 0.10 |b/in.? are given in table 3 for information and 
eorrespond to blankets having different compressibility. 





Properties of Household Blankets 283 


TaBLE 3.—Minimum thickness at 0.10 lb/in.? suggested for blankets of various 
weights and compressibilities 








in./in 
Compressibility, 
Ib/in 












































TaBLE 4.—Thermal transmission and thermal insulating values corresponding to 
various thicknesses at 0.10 lb/in? 








Thermal 


Thermal 
: insulating 
transmission value 








Btu/(°F 
hr ft*) 
1, 52 





45 
. 33 
. 23 
-15 
. 08 
- 00 
. 93 
. 87 
81 
yf 
- 66 














5. BREAKING STRENGTH 


The breaking strength was found to be related to the weight and 
to the compressibility of a blanket. Minimum values for breaking 
strength (grab method) of the warp and filling for blankets of various 
weights and compressibilities are suggested in table 5. Manufac- 
turers should have no difficulty in meeting the strength requirements 
suggested. Most blankets meeting the other requirements will ex- 
ceed the suggested values for breaking strength. Furthermore, since 
blankets generally increased in breaking strength during laundering, 
it is believed that the values suggested in table 5 will assure satisfac- 
tory performance of the blankets for ordinary use. 


586263—44———5 
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TaBiEl 5.—Minimum breaking strength (grab method) suggested for blankets of 
various weights and compressibilities 








Warp strength Filling strength 





./in. 
/in.2 


j . Pasa Pe, | 
Weight Compressibility, ibjin | Compressibility, Ib 





Below | 0.25to | Above | Below | 0.25to | Above 
| 0.35 | 0.35 0.25 0.35 | 0.35 

















RSSSSSERSS 














6. AIR PERMEABILITY 


No requirements are suggested for air permeability, since in general 
blankets will be used in nominally still air. For outdoor uses 
a combination of blanket-and-wind resistant cloth has decided 


advantages.’° 
7. SHRINKAGE 


The maximum allowable shrinkage suggested is 5 percent in the 
warp and in the filling. This amount of shrinkage is to be expected 
in normally manufactured blankets. Proper precautions in laundering 
blankets to preserve their characteristics should automatically 
safeguard against excessive and objectionable shrinkage. It seems 
advisable to specify nominal dimensions of the blankets 5 percent 
greater than necessary so that the blankets will not be too small 
after a shrinkage of 5 percent. 


Wasuineton, December 28, 1943. 


10 Herbert F. Schiefer, Adrantages of a blanket-and-sheet combination for outdoor use, J. Research NBS 
30, 209 (1934) RP1529. 
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STUDIES OF PORTIONS OF THE QUATERNARY SYSTEM 
SODA-LIME-SILICA-WATER AT 25° C 
By George L. Kalousek 


ABSTRACT 


A study has been made of portions of the system soda-lime-silica-water at 25° C, 
The only solid bor ere found were Ca(OH), and a four-component gel of variable 
composition. he boundary, Ca(OH),—soda-lime-silicate gels, was determined 
and the compositions of the gels along this boundary shown to vary from 
0.003Na,0 :2.0CaO:1.0Si0O,:2H,O (at 0.2 g of Na,O per liter) to about 0.25Na,0:- 
1.0Ca0:1.08i0,:2H,0 (at 20 g per liter and extending to 101 g per liter of Na,O). 
In regions off the boundary, at selected constant concentrations of Na,O but with 
increasing concentrations of SiO, in solution, the Na,O:SiO, molar ratio of the 
gels varied only slightly from 0.2 in most cases; the CaO:SiO, molar ratio, 
however, decreased to values approaching 0.1 at the maximum concentrations 
of SiO, used. Interpretations pertaining to relations between the composition of 
the gels and solutions are given. 


CONTENTS 
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. Materials 
. Procedure 
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2. Methods of chemical analyses 
. Results and discussion 
1. Preliminary observations 
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(a) Compositions of soda-lime-silicate gels__.......____- 289 
(b) Relation between compositions of gels and solutions 
along boundary 
3. Regions of higher silica concentrations 
. Applications of the results of the study to concrete 
. Summary 
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I. INTRODUCTION 


The role of the alkalies, soda and potash, in hydrating cements is 
not fully understood. Information on how these constituents alter 
the chemical composition of the lime silicate hydrates and of other 
products in hydrating cements should be highly useful in explaining 
the mechanism of the various hydration processes involved. Also, 
a knowledge of the compositions of the alkali solutions in which the 
hydrated cement compounds are stable should prove helpful in eluci- 
dating the effects of the alkalies on hydrating cements in the presence 
of “reactive” aggregates. 


285 








286 Journal of Research of the National Bureau of Standards 


The present investigation was accordingly designed for obtaining 
a portion of such data and was concerned with the effects of Na,O on 
the compositions of the lime silicate hydrates and of the aqueous phase 
in which they are stable. The investigation was rth a study of 
portions of the quaternary system soda-lime-silica-water at 25° C. 


II. MATERIALS 


The starting materials were solutions of sodium silicate and sodium 
hydroxide of known compositions, crystals of Ca(OH), and distilled 
water. 

The sodium silicate solution was a reagent-quality product con- 
taining not more than 0.2 percent of CO, and having a molar ratio of 
1.00Na,0:3.34SiO,. Sodium hydroxide was freed of carbonate by 
preparing a saturated solution. The clear supernatant liquid was 
drawn off and diluted as required. 

Calcium hydroxide of reagent quality and having a particle size of 
less than 1 micron was used for most of the experiments. Also, 
Ca(OH), crystals of varying sizes, prepared by adding solutions of 
NaOH to saturated solutions of Ca(OH)2, were used in a number of 
exploratory preparations. 


III. PROCEDURE 
1. PREPARATION AND FILTRATION OF MIXTURES 


A few tentative mixtures of the standard solutions and the 
Ca(OH), were made in order to establish the amounts required to 
produce the desired quantities (about 0.5 to 1.0 g) of the precipitates 
in mixtures having a volume of 200 ml. Most of the mixtures were 
made by adding a solution of sodium silicate to a well-agitated sus- 
pension of fine-grained Ca(OH), in a solution of NaOH. The result- 
ing preparations in tightly stoppered ceresin-lined flasks were stored 
in an air bath maintained at 25+0.05° C until filtered at ages rang- 
ing from 2 to 18 weeks. During the first 5 to 8 hours after mixing 
the preparations, they were shaken about once each hour, and then 
once daily after the first day. Although this was the general proced- 
ure followed, it will be necessary, as required in certain sections of 
this report, to give more specific details. 

The time allowed the mixtures to ‘age’ was selected arbitrarily. 
It may be mentioned, however, that neither the solid phase nor solu- 
tion underwent any determinable changes in time periods of from 
several weeks to several months. It is recognized that this steady 
state may be displaced by eventual crystalization of the solid phase, 
a process which, however, appears to proceed exceedingly slowly at 
room temperature. 

The mixtures were filtered through a medium-textured sintered glass 
crucible by application of suction. The residue on the filter during the 
last stages of filtration was firmly compacted with a rubber plunger, 
and as much of the solution was withdrawn as permitted by the 
retentivity of the gelatinous precipitate. The precipitate was next 
redispersed in a small amount of wash solution (ethyl alcohol and 
water or water if the filtrate contained more than 70 g of SiO, per 
liter) and the solution again filtered off. The most suitable ratio of 
alcohol to water for each wash solution was determined by testing a 
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few milliliters of the original filtrate and selecting the maximum alcohol 
content which just failed to produce a suspension. Four to five 
successive washings were made with solutions of progressively increas- 
ing alcohol content, then four to five washings with 95-percent ethyl 
alcohol, and finally four to five washings with ethyl ether. In the 
course of each washing, the precipitate was redispersed and during 
filtration, compacted on the filter. After being washed the precipitate 
was left in the laboratory air for 1 day and then stored in a tightly 
stoppered vial. 


2. METHODS OF CHEMICAL ANALYSES 


Standard procedures of chemical analyses were followed. All of 
the solutions and most of the precipitates were analyzed for soda, lime, 
and silica. The “hydrate” water of the precipitates was determined 
by difference. 

The amounts of Na,O in the filtrates containing relatively small 
amounts of SiO, (1 percent or less of the total solids) were determined 
by titrating the solutions with HCl, using methyl red as indicator. 
Allowance was made for the Ca(OH), present, which was determined 
gravimetrically. The errors in the results were not greater than 1 
percent of the amount of Na,O, as determined occasionally for check 
purposes by precipitation. In all the other determinations, the 
Na,O, after removal of SiO., was precipitated and weighed as sodium 
zinc uranyl acetate hydrate. 

Calcium was precipitated as oxalate, ignited, and weighed as CaO. 
However, when the filtrates contained over 50 g of Na,O per liter the 
CaO was precipitated as phosphate, converted to sulfate, and reprecip- 
itated as oxalate according to the method described by Hillebrand 
Lundell [1].! In the analysis for SiO, the solutions were dehy- 
drated twice with HCl. The combined residues were igmted to con- 
—_ weight and treated with HF. The loss in weight was reported 
as Si 2° 

IV. RESULTS AND DISCUSSION 


1, PRELIMINARY OBSERVATIONS 


Exploratory mixtures showed that only gelatinous precipitates were 
formed. Some of these, in mixtures of high concentrations of SiO,, 
were highly solvated gels and others, in mixtures of high concentrations 
of Na,O and low concentrations of SiO., appeared as grains of gel. 
Analyses indicated that the CaO:SiO, molar ratio of the gels decreased 
as the concentration of Na,O in solution was increased and the 
Na,O:SiO, molar ratios were found to increase up to a certain value. 
At constant concentrations of Na,O in solution, but with increasing 
amounts of SiO., the CaO:SiO, ratios were found to decrease. The 
four-component gelatinous products are referred to as soda-lime-silicate 
gels, or, more often, simply as gels. 

The shifts in the compositions of the gels could not be followed 
microscopically as the changes in the indices were small and appeared 
to be affected by aging. The presence or absence of small amounts of 
fine-grained crystals of Ca(OH), could not be ascertained reliably by 
means of the microscope. Large crystals of Ca(OH), could not be 


Figures in brackets indicate the literature references at the end of this paper. 





288 Journal of Research of the National Bureau of Standards 


used as they were converted but slowly to the siliceous products owing 
to the formation of gelatinous deposits around the crystals. 


2. BOUNDARY, Ca(OH):,—SODA-LIME-SILICATE GELS 


Considerations based on the preliminary results indicated that the 
best initial procedure to follow in studying the system would be to 
ascertain the compositions of the aqueous phase in contact with 
Ca(OH), and the soda-lime-silicate gels. Accordingly, the mixtures for 
establishing this boundary were prepared, only a moderate excess of 
crystalline Ca(OH), being used. The concentration of Na,O was varied 
from 0.2 to 152 g per liter. After 2 to 7 weeks, the mixtures were 
filtered and the filtrates and precipitates analyzed. 

All the data pertaining to the mixtures, filtrates, and precipitates 
are reported in table 1, the compositions of the solutions being given 
in grams per liter of Na,O, CaO, and SiO,. The table shows that the 
amount of CaO in solution decreased from 0.92 to about 0.02 g per 
liter as the concentration of Na,O was increased from 0.2 to about 
20 g per liter. Simultaneously, the concentrations of SiO, increased 
from 0.002 to about 0.020 g per liter. In the range of concentrations 
of Na,O from 20 to 152 g per liter, the amounts of CaO in solution 
remained nearly constant at a value of about 0.010 g per liter, but 
those of SiO, increased from 0.02 to over 1.1 g per liter. 

Figure 1 shows the g per liter of CaO and SiO, found in solution at 
the different concentrations of Na,O. 
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Ficure 1.—Amounts of Na,O, CaO, and SiO, in filtrates from mixtures containing 
solid Ca(OH), and soda-lime-silicate gels. 
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TaBLE 1.—Compositions of solutions in contact with crystalline Ca(OH), and soda- 
lime-silicate gels and the compositions of the mized solids 



































Composition of Mix- Composition of solu- : 
tures in grams per liter . tions in grams per liter| Molar compositions of mixed solids 
ge 
No in | g 
— weeks | From direct analyses 
Na20*| CaO | SiO: Na:0 |} CaO | SiOs | Computed | 
CaO: SiOs l 
Na20; CaO | SiO:| H:0 
one 2.0 | 1.864 | 0.655 4/ 1.69 | 0.380 | 0.001 2.43 | 0.03 | 2.42 | 1.00 j_..... 
, oa 3.5 | 1.490 655 4 3.16 - 195 . 001 2.12} .05 | 208 | 1.00 |...... 
RR ee i 5 4.5 | 1.248 . 655 7 4.00 . 140 . 006 b1.83 -05 }-1.00 | 1.00 |...... 
Pa Dente 6.5 | 1.116 - 655 4 6.17 . 090 . 002 1. 69 sini het 1,00 t..... 
§.. woah yg 0, 932 . 655 4 8.74 055 . 002 1, 43 -13 | 1.38 | 1.00 3.3 
PG a rts 13 . 745 . 655 4 | 12.00 037 .O11 1.18 -18 | 1.18 | 1.00 3.0 
ee ‘ 16 . 947 . 655 7 | 15.30 029 . 008 1.52 +e eee Fes | 
Ye ee 18 . 655 640 4/| 17.8 019 015 1.09 . 22 | 1,09 | 1.00 3.3 
i cede 21 . 618 640 4 | 20.6 .013 - 026 1.06 22 | 1.03.) 1.00 29 
10 2 | .806| .655| 4/243 | :009| .043 1.03 | .23| 102/100} 3.1 
| SRR RT eat 37 1.490 | 1.310 6 | 36.5 . 038 . 046 1.23 23 1.25 | 1.00 2.9 
12 a eee 40 0.945 | 0.655 2} 39.0 . 028 . 048 b1.63 | .29) 1.31 | 1.00 3.6 
Zeon 40 - 851 655 5 | 39.6 025 . 065 RP de dand Hincnattcoccskscomp 
14... aes 55 1.490 | 1.310 6 | 54.5 018 . 071 1. 24 26 | 1.26 | 1.00 3.2 
15 1.490 | 1.310 6 | 67.0 oll . 106 1,32 25 | 1.34 | 1.00 3.0 
) RRO BoE 77 0.605 | 0.655 2 | 76.2 008 . 325 boel 23 | 1.08 { 1.00 3.0 
| RE: 78 - 606 5 | 77.4 004 - 227 i Seen Tha FRG Sy Fae 
i ecanasnt 82 1.490 | 1.310 6 | 80.8 -010 . 157 1. 38 .25 | 1.37 | 1.00 3.0 
Wl kscebenss 100 1,490 | 1.310 6 | 99.0 - 004 . 241 61.49 .24 | 1.36 | 1.00 2.9 
Oe dices 110 0,932 | 1.310 7 |108.1 . 006 .377 1.07 27 | 1.02 | 1.00 2.7 
| EMR Seep. 110 . 756 | 1,192 2 /109 . 004 - 428 1,05 24 | 1.06 | 1.00 3.0 
| Se es 116 - 932 | 1.310 6 |115.0 . 008 . 424 ee eee 1,12 | 1.00 2.7 
, aan eels 125 - 932 | 1.310 6 |123.5 .010 . 501 12 .29 | 1.26 | 1.00 3.2 
> SER eee 135 - 932 | 1.310 6 |132.5 . 006 . 633 1. 46 .25 | 1.43 | 1.00 3.2 
Geiniaccaenen 144 . 880 | 1.310 3 (141.6 .005 | 1.058 1. 60 t Bet SE" g Ree 
26 --| 145 . 932 | 1.638 5 |143.5 .006 | 0.845 1,27 . 26 | 1.26 | 1.00 2.9 
| SEER Re et> . 932 | 1.965 5 152.0 .007 | 1.131 1.18 - 2 | 1.16 | 1.00 2.6 
Wonka mecawatel 155 . 530 | 1.310 5 |153.1 .00 Sie Ureccene<seue A Sage Ig ES ARE See 4 
| 





























® The amount of Na;O added is approximate in these preparations. 
> Ca(OH): separated partly from mixture of solids. 


(a) COMPOSITIONS OF £ODA-LIME-SILICATE GELS 


The precipitates obtained in the mixtures used for locating the 
boundary, although containing Ca(OH)s, were suitable for ascertaining 
the Na,O:SiO, molar ratios of the soda-lime-silicate gels. 

The compositions of the precipitates, assigning a value of unity to 
the mols of SiO, are given in table 1 (last four columns) as the molar 
ratios of NagO:CaO:SiO,.:H,O. Thus, the column under Na,O gives 
the molar ratio of Na,O:SiO;, of the gels. 

The molar ratios of CaO to SiO, of gels along the boundary could 
not be determined reliably by preparing mixtures, computed to be 
removed just slightly off the boundary so as to avoid the presence of 
Ca(OH),, and analyzing the precipitates. Because of this difficulty, 
the following procedure for determining the composition of the gel 
coexisting with crystalline Ca(OH), at constant concentration of 
Na,O was selected as the most expedient and the one likely to give 
reliable results. 

Series of 5 to 10 mixtures each were prepared at Na,O concentra- 
tions (or levels) of 0.2, 2.4, 5.5, 9.2, 11.6, 22.3, 51, 75, and 101 g per 
liter. (See table 2.) The total quantities of Ca(OH), and SiOz, 
added as sodium silicate, were estimated from the amounts needed 
for the solutions as shown in figure 1, plus the approximate amounts 
required for the precipitates as indicated by results on exploratory 
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TaBLE 2.—Data 
coexisting with 


rtaining to the CaO: SiO, molar ratios of soda-lime-silicate gels 
a(OH), and solutions containing between 0.19 and 101 g of Na,0 




































































per liter 
Composition of mix- Composition of solu- 
tures in grams per tions in grams per | Composition y)> »-ppmmaane molar 
liter liter 
Age, — 
No. es 7 oe Direct analyses 
lompu 
Na:O0 CaO SiO: Na20 CaO SiO, CaO: SiOs | 
Na0| CaO 810s | 10 
1.....-.--.--| 0.205 | 1.891 | 0.655 6 | 0.195 | 0.862 | 0.002 hf aaeed Ae a 
Nao .200 | 1.968} .655]| 18| .190| .888| .002 ky Re apres cape se 
Rye 205 | 1.968 | . 655 6| .195| .858| .002 Sian) ...5 meth 
4.......---.-| .205 | 2.042 | .655 6| .195| .886| .002 tf Ree See Rs, eens i: 
sk ever .205 | 2.119} .655 6| .195| .912| .003 1.98 |0.003 | 1.94 | 1.00 | 3.2 
CORRELA .205 | 2.192 | . 655 6} .190| .894] .002 hy harman set 
_ BE Shite .200| 2.270} .655|} 18] .186| .915| .002 Se Raa Wigs tothe 
WRIA .200| 2.420] .655] 18] .188| .911| .002 ST. huss c. ce eae 
9.......--...| 92.5 | 1.058] .655 5| 2.37 | .193| .003 oy Rea coe | 
RC. 2.5 | 1.135 | .655 5] 2.38 | .237| .003 1.47] .04 | 1.41 | 1.00} 33 
RR eat 2.5 | 1.210] .655 5| 2.36 | .240| .003 Lae ea 
1$......... 1 2S 12a ee 5| 2.38 | .243| .003 Yh eR aE EINE cata 
RISER T 2.5 | 1.361] .655 5| 238 | .245| .002 CA SM ear sadam, LI Oe aN 
a Se ts) 5.8 | 0.525 | .655 5| 5.40 | .002| .039 0.92 | .18 | 0.90 | 1.00 | 2.60 
AI 5.8 600 | . 655 5| 5.42 | .010] .015 .99 1.15 | .98| 1.00] 268 
wate GP 675 | .655 5| 5.42 | .022| .009 1,08 | .13 | 1.07 | 1.00 | 2.80 
Sp Pa, 5.8 750 | . 655 5| 5.42 | .041| .005 1.17] .12 | 1.15 | 1.00 | 3.00 
a al 5.8 826 | .655 5| 5.46 | .063| .005 1.24] .11 | 1.21 | 1.00 | 3.06 
RES Beals 8 5.8 .900 | .655 5| 5.46 | .100] .005 1.32 | .07 | 1.30] 1.00 | 2.98 
at pce. 5.8 975 | .655 5| 5.46 | .099| .003 1.44] .08 | 1.41 | 1.00 |(3.10) 
ARB a 5.8 | 1.050] .655 5| 5.46 | .101| .002 1.55 | .08 | 1.52 | 1.00 |(3.35) 
ES i 5.8 | 1.125] .685 5| 5.46 | .007| .003 4 ORES Gita ako hale 
aR EA aS 5.8 | 1.200] .655 5| 5.52 | .103| .003 ¥, ) RDN sealpes neal 
ees 9.5 | 0.522] .655 6| 9.14 | .011| .060 UN Ne Se SERS Sa aa 
ee 596 | . 655 6| 9.15 | .015| .038 i> ADIN RR IBS, 5 2 OR 
RNP ee 9.5 | .671| .655| 6| 9.17 | 1025| ‘017 ey Fi ae peter Lee 
a 745 | 655 6| 9.15 | .039| .012 1.18 | .14 | 1.13 | 1.00 | 29 
eae 820 | .655 6| 914 | .046| .007 | Ng PCG eile Sete 
ieee 8 9.5 .894 |] .655 6| 9.14 | .044 as LOD dissickcdswisnsiessuc 
_ eee 530 | .655 5| 11.55 | .003| .060 OU elie siet Baath fen. voit 
RIN 2 12.0 606 | . 655 5| 11.53 | .007| .025 Fy a Bilctie shatbin: bite 
ahaa? 12.0 682 | .655 5 | 11.56 | .023| .014 1.05 | .19 | 1.05 | 1.00 | 2.7 
eee 757 |. 655 5 | 11.56 | .036| .008 J) 3, LOS ee treats 
Diiasiainstite 12.0 833 | .655 5 | 11.60 | .036| .008 5 DSS ara Vialagee - 
RS cae 22.0 .873 | .655 5 | 22.11 | .002] .236 Urn Mae 
36...........| 22.0 1466 | | 655 5 | 22.11 | :003| .149 ecches- 
mb 559 | .655 5 | 22.30 | .016] .071 1.00 | 2.6 
i...) ae 652 | .655 ST Tt el at” CO a od ‘wan 
PRB, FY 746 | 655 $49008' 401i 0004)2 © 1.80 telichkosdacal 
40....../....| 51.0 | 1.092/1.705| 8/508 | .o1 | .488/ 0.05|-_..|...-. J}. 
ya Bes 51.0 | 1.390 | 1.705 8 | 50.9 /01 | .178 1.00 | 2.7 
RES Bias, 51.0 | 1.688 | 1.705 8 | 50.7 /01 | .105 silks & 
43.......----| 51.0 | 1.985 | 1.705 8 | 50.5 fats eee. (240 i oles oes aie 
eR 51.0 | 2.273 | 1.705 8 | 50.5 /01 | .070 -+| : 
~ ae ae 75.2 | 0.993 | 1.772 4| 75.4 a Agee Sd 
a PEM “| 75.2 | 1.242} 1.772 rH haalec Ga we ihe ‘on 8 
PERI 75.2 | 1.490 | 1.772 4 | 75.2 01 | .220 | 1.00 | 2.9 
48...........| 75.2 | 1.787|1.772| 4| 75.3 | ior | :121 | ts eae 
Boidctias 75.2 | 1.985 | 1.772 4 | 75.5 01 | .106} ° AY onset 
101.5 | 0.993 | 1.772 4 |101.2 .01 | .721 hk Ss 
“/101.5 | 1.242 | 1.772 4 |100.6 01 | .498 1.00 | 2.6 
~.|101.5 | 1.490 | 1.772 4 |100.8 01 | .320 | oe 
--{101.5 | 1.737 | 1.772 4 |101.0 01 | .289 | 
~.---|101.5 | 1.985 | 1.772 4 |100.8 .01 | .350 | ie 
* The amount of NaO; added is approximate in preparations 9 to 54. 
mixtures. The CaO:SiO, ratio was varied systematically in each 


series so that, after reaction, about half of the mixtures contained 


Ca(OH), and the remainder none. 


After 4 to 18 weeks the mixtures 


were filtered and the filtrates analyzed. The computed Ca0O:Si0; 
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molar ratios,’ representing the differences between the amounts of 
CaO and SiO, in the mixtures and the amounts found in solution, 
are used for the graphical presentation of the data. Most of the 
precipitates in the series containing 5.4 to 5.5 g of Na,O per liter were 
analyzed for check purposes, and in each series those having composi- 
tions nearest to the boundary curve were also analyzed. The re- 
sults are given in table 2. 

The CaQ:SiO, molar ratios of the precipitates at the different Na,O 
levels are plotted in figure 2 against the amounts of CaO or SiO, 
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Ficure 2.—Plots of the CaO:SiO, molar ratios of precipitates against the amounts of 
CaO or SiO, in solution for each of nine series of mixtures of varying CaO and SiO, 
contents but of constant Na,O concentrations. 

O—O=Ca0. 


@....@=Si0;. 


(or both) in grams per liter of ‘solution. The CaO:SiO, molar ratios 
of the gels coexisting with Ca(OH), at different Na,O levels are those 
corresponding to the “‘breaks”’ in the curves shown in figure 2. These 
breaks were quite sharp except for the series of mixtures containing 
small amounts of CaO or SiO, in solution. This may be due to the 
experimental errors in the determination of these small quantities and 
therefore, to obtain as nearly a representative value as possible, the 
respective concentrations of both constituents were plotted. The 
average CaQ:SiO, ratio of the two was selected as being more accurate 
than either one. 
The approximate water content of the gels, as moles of H,O per mole 
of SiO,, is given in table 2. However, in table 1 the water content 
? Both the computed and observed CaO:SiO: molar ratios were determined, and are given in table 1 for 
the purpose of indicating the reliability of the former. With a few notable exceptions, in which the ob- 
served ratios were markedly lower than the computed ratios, the two agreed within a few hundredths of 
a mole. Such small differences are not serious, and knowing the cause for the exceptions, the computed 
ratios | were used in some of the following parts of the study. The marked exceptions, which occurred 
in the earlier experiments, were due to a partial separation of the solid phases, the Ca(OH)s adhering firmly 
to the ceresin ous of the flasks. In subsequent mixtures such separation was prevented by using a finel 


crystalline Ca(OH)s and vigorously shaking the mixture during preparation and also at frequent interva 
during the first day. The resulting suspensions did not cake in the bottom of the flasks. 
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of the precipitates showing a CaO:SiO, molar close to 1 (the precipi- 
tates being mostly gel) also indicates approximately the water con- 
tent of the gels. 

The molar ratios of Na,O0:SiO,, Ca0:SiO,, and H,O:SiO,, repre- 
senting the compositions of the gels along the boundary, Ca(OH), — 
soda-lime-silicate gels, are presented graphically in figure 3. This 
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FicurE 3.—Relations between the concentrations of Na,O in solution and the com- 
positions of corresponding soda-lime-silicate gels expressed in molar ratios of 
H,0:SiO,, CaO:SiO,, and Na,O:SiQ,. 


figure, based on data from tables 1 and 2, and figure 2, shows that the 
Na,O:SiO, molar ratio increased from a value of 0.003 to a-value of 
about 0.23 as the concentrations of Na,O in solution were increased 
from 0.2 to about 20 g per liter and them remained nearly constant 
* this value for all higher concentrations of Na,O, up to 152 g per 
iter. 

The CaO:SiO, molar ratios of the gels, as shown by figure 2, de- 
creased from 1.98 to 1.00 as the concentration of Na,O was increased 
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from 0.2 to 22.3 g per liter. At the higher concentrations of Na,O 
used the ratio remained nearly constant at a value of about 1.0. 

The H,0:SiO, molar ratios varied between 2.5 and 3.3. This 
spread in the results is rather large, but not surprising, considering 
the nature of the materials. 


(b) RELATION BETWEEN COMPOSITIONS OF GELS AND SOLUTIONS ALONG 
BOUNDARY 


In order to show the relation between the compositions of the gels 
and solutions along the boundary, the data for mixtures 5, 10, 19, 27, 
32, and 37 in table 2 were computed as mole percent of Na,O, CaO, 
and SiO,, the water being disregarded. These results are plotted in 


5/02 
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MOLE PERCENT 
FIGURE 4.—Compositions of gels coexisting with Ca(OH), and of the contact solutions. 
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The results are expressed in mole percent of Na3O, CaO, and SiO; (anhydrous basis), 


O=Gel. 
@=Solution. 


THEORETICAL 


@=).25Na2,0.Ca0.SiOx 
+=—=2Ca0.SiOx. 


figure 4. It should be noted that gels and solutions obtained at Na.O 
concentrations above 22 g per liter had an average composition, in 
mole percent, almost identical with that of mixture 37 table 2, and 
therefore would fall nearly on the same point. 

Figure 4 shows that the amounts of SiO, in solution along the 
boundary were relatively small throughout the entire range. It is 
seen, however, that the decreases in the amounts of CaO in solution 
caused large changes in the compositions of the gels. But, once the 
CaO concentrations attained an approximately constant value, the 
compositions of the gels remained nearly constant. 
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3. REGIONS OF HIGHER SILICA CONCENTRATIONS 


The effect of the SiO, contents of the solutions on the compositions 
of the gels was studied by maintaining the concentrations of Na,O at 
approximately 6, 13, 50, and 96 g per liter, and increasing the concen- 
trations of Sid, in each series to a value of 3.3 moles of SiO, per mole of 
Na,O. These mixtures contained insufficient Ca(OH), to saturate the 
solutions. 

Both the solutions and solids were analyzed and the results are given 
in table 3. The compositions of the gels were computed only for the 
6- and 13-g Na,O levels, since the computations for the 50- and 96-g 
Na,O levels would be subject to large errors due to the relatively 
small amounts of constituents that precipitated compared to the 
amounts that remained in solution. The occasional large difference 
between computed and observed compositions at the 6- and 13-gram 
Na,O levels may be due to the same cause. 

The last nine columns of table 3 give the compositions of the gels 
and solutions in mole percentages of Na,O, CaO, and SiO,. It may be 
mentioned that the compositions of the gels computed to mole percent 
of the four constituents showed that the H,O content was nearl& 
constant at a value of 55 +3 mole percent for all but two of the gels at 
the different Na,O levels. It is evident, therefore, that presenting 
the data for the gels as a ternary system is in reality showing the 55 
mole percent H,O-level of the quaternary system. 

The results for the 6-g and 13-g Na,O levels are plotted in figure 5, 
and those for the 50- and 96-g Na,O levels in figure 6. The effect of 
Na,O on the compositions of the solids is readily ascertained by com- 
paring the diagrams in the two figures. In order to facilitate the 
comparison, the amounts of SiO, in solution, in grams per liter, are 
given in the diagrams to the right of the points representing the com- 
positions of the solutions. 

Both the computed and observed compositions of the solids are 
plotted in figure 5. Considering the highly solvated nature of the 
gels, which made filtration ro | washing difficult, and that small 
analytical errors caused large errors in the computed results, the 
differences between the two sets of values are not to be regarded as 
significant. The trends of both, however, are about the same, and 
the averages may be taken as representative of the true compositions 
of the gels. Figure 5 shows that as the SiO, in solution was increased, 
the mole percent of Na,O in the gel did not vary greatly, but that of CaO 
decreased markedly. A comparison of the two diagrams in figure 5 
indicates that the gels in contact with solutions of the same mole 
percent composition have about the same compositions. Noting that 
a given mole percent of SiO, at the two Na,O levels does not represent 
the same concentration of SiO, in solution, it follows, therefore, that 
the ratio of SiO,:Na,O (or Na,O:SiO,), in solution governs greatly the 
compositions of the gels. 

About the same remarks apply to figure 6 as were made for figure 5 
except that the curve for the 50-g level, above 65 mole percent of 
SiO, in solution, is displaced slightly further to the right than those 
at the 6- and 13-g¢ Na, O levels. ether this is an actual shift in 
compositions or is due to an incomplete washing out of the Na, 0 
and SiO, in the solution associated with the gel is not known. The 
mixtures at the 96-g Na, O level, having SiO, concentrations in excess 
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of 90 g per liter, contained exceedingly fine-grained precipitates which 
passed readily through the filter even after 7 weeks of aging and, 
therefore, only a few filtrates were available for analysis. The few 
results available at this level do not show as smooth a relation as those 
at the other levels. 


Si0p 
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Figure 5.—Compositions of two series of solutions containing 6 and 13 g of Na,O 


per liter, respectivel,y and of the corresponding gels in mole percent of NazO, CaO 
and SiO, (anhydrous basis). 





O=Gel, analyzed, 
X = Gel, computed. 
@= Solution, analyzed, 
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V. APPLICATIONS OF THE RESULTS OF THE STUDY TO 
CONCRETE 


This investigation has shown that gels coexisting with crystalline 
Ca(OH), and with solutions containing between 0.0 and 100 g of 
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FicuRE 6.—Compositions of two series of solutions containing 50 and 96 g of Na,;0 
per liter, respectively, and of the corresponding gels in mole percent of Na,O, Ca0, 
and SiO, (anhydrous basis). sai 

=Solution. 


Na,O per liter have compositions that range between about 
2CaO.1Si0,.czH,O and 0.25Na,0.1Ca0.1Si0O,.2zH,O. If similar rela- 
tions prevail in a completely hydrated cement, the approximate com- 
positions of the gel and solution may be deduced as follows: 
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For purposes of illustration, a cement containing 20 percent of 


” SiO, and 1 percent of Na,O is used. The extent of combination of 

Na,O and SiO, as a soda-lime-silicate gel will depend on the extent 
ne to which the Na,O-bearing compounds and the two lime silicates have 
of hydrolyzed. Because of this the simplest deductions, as follows, are 


those for a completely hydrated cement. If all of the Na,O combined 
with SiO, as soda-lime-silicate gel, the Na,O:SiO, molar ratio of the 
gel would be 0.05. As shown in figure 3, a gel of this Na,O:SiO, ratio 
is stable in contact with a solution containing 3 g of Na,O per liter. 
Since this amount * of Na,O in solution is only 2 percent of the total, 
it may be disregarded in calculating the composition of the gel. 
Figure 3 also shows that a soda-lime-silicate gel coexisting with a 
solution containing 3 g of Na,O per liter has a CaO:SiO, molar ratio 
of 1.4. A combination of the two ratios gives the following composi- 
tion for the gel, 0.05Na,0:1.4CaO:1SiO,:2H,O. The solution in 
contact with this gel would contain about 3 g of Na,O, 0.1 g of 
CaO, and 0.005 g of Si0, per liter. 

It is recognized that the above calculations are not entirely valid. 
For example, it is probable that a portion of the silica in a completely 
hydrated cement is present in combinations other than soda-lime- 
silicate gel. Also the composition of the gel may be altered by the 
prolonged aging required to bring about complete hydration of the 
cement. 

In a partially hydrated cement, the compositions of the gel and 
solution will differ from those given above. It is known that the 
dicalcium silicate in a cement paste is only slightly hydrated at 28 
days, but the tricalcium silicate in the same time may approach a 
state of complete conversion to its hydration products. Since, on the 
average, about one-third of the SiO, in a cement is combined as 
dicalcium silicate, it would follow that, the gels formed during the 
earlier periods would contain more Na,O than those prevailing at 
complete hydration, providing the Na,O-bearing compounds react 
with water at rates comparable to that of tricalcium silicate. Also the 
solution would be richer in soda. Making allowance for the average 
value of SiO, combined as dicalcium silicate, calculations indicate 
that the gels formed during the earlier periods may have Na,O:SiO, 
molar ratios as high as 0.1, and the contact solution may contain as 
much as 10 g of Na,O per liter. 

The different values reported in the literature for the compositions 
of the silicate gels in hydrated cements may be related to the alkali 
contents of the cements as indigated by the following considerations. 
Bogue and Lerch [2], in a study of the alkali-free lime silicates, came 
to the conclusion, from results of free-lime determinations, that the 
hydrated products had compositions approaching a value of 2CaQO: 
$i0,:2H,O. This result is in accord with that obtained in studies of 
the system line-silica-water reported by Flint and Wells [3] and 
Bessey [4]. Using cements that contained between 0.6 and 1.5 
percent of total alkalies, Bessey [5] concluded from calorimetric 
results and free-lime determinations that the hydrous silicates formed 
in the pastes had a composition of about 3CaO:2SiO,:zH,O. This 
value of the CaO:SiO, ratio is in fair agreement with that computed 

+ In making these computations, the C/W ratio was selected as 2.86 by weight, and it was assumed that the 


“fixed”? water amounted to 30 percent by weight of the cement. No allowance was made for water gained 
or lost by the paste. f 
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~ - study for hydrated cements containing about 1 percent of 
a0. 

Concretes made from high-alkali cements and aggregates containing 
opal or certain other active siliceous materials have been shown b 
numerous investigations to be unsound. A gelatinous product formed 
in such concrete is generally associated with this unsound state, 
Analyses of some of these gels show them (see table 4) to consist 
largely of silica, alkalies, and water (ignition loss), together with 





V/ V V V \V/ V 
CaO MOLE PERCENT 


Fiaure 7.—Compositions of synthetic gels and of gels found in concretes showing 
expansion and cracking 


The results are expressed in mole percent of NazO, CaO, and SiO» (anhydrous basis). 


@—@ Composition of gels, solution containing 6 g of Na2O/liter. 

x——-xX Composition of gels, solution containing 13 g of NasO/liter. 

O——©O Composition of gels, solution containing 50 g of Na2O/liter. 

A—~—A Composition of gels, solution containing 96 g of Na2O/liter. 
--- Composition of gels coexisting with Ca (OH): 





COMPOSITION OF GELS FROM CONCRETE 


= Meissner. C= Meissner. 

= Coombs. X=Blanks. 
A=Stanton. = Meissner. 
V =Stanton. = Hanna. 


*= Average of above. 


smaller quantities of CaO, R,O;, MgO, and SO;. The data computed 
to mole percent of Na,O (the K,O being counted as Na,O), CaO, and 
SiO, are presented in figure 7. Since some of the constituents were 
disregarded in these computations, the results are to be considered as 
only approximating the true ratio of Na,O to CaO to SiO, of the 
gels. Giving consideration to this and the fact that the gels were 
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obtained from different sources under different circumstances, the 
spread in the data is not surprising. The average of these composi- 
tions is given in figure 7 for facilitating comparisons between the 
synthetic gels and those found in concrete. The compositions of the 
gels shown in figure 7 are those previously given in figures 4, 5, and 6. 


TABLE 4.—Percentage oxide compositions of gels found in or on concretes contain- 
ing ‘‘reactive’’ aggregates 














Number | SiO» R203; | CaO | MgO SO; | Na,O | K;0 Loss | Reported by— 

he 81.9 0.7 Ey a BS a 4.0 2.6| 10.5 | Meissner [6] 

Qa. 85.0 2.2 eS Eten ers | 4.4 1.9 16.1 | Coombs [7 

3 53. 4 7.8 2.6 eh eee | 128 b ¢16.7 | Stanton [8 

44 53.9 6.8 2.9 AY Sea 12.9 i © 22.4 Do. 

a; 53. 4 1.5 4.7 .0 0 12.7 4.5 20.5 | Meissner [6] 

€.5.. 58.7 5 0 0 0 19.9 6.5 13.8 | Blanks [9 

fe EA Linwancth Gngiaeion sin calibre 0 16.8 ae Ee Meissner [6] 
Re Sas 50.0 5 2.8 2 2.6 24.7 8.2 11.0 | Hanna [9] 





























« Analysis of water soluble portion of gel amounting to 90.9 percent of total (ignited basis)— loss on igni- 
tion determined on gel as found. 

b K30 counted with Na,O. 

e Sum of “Moisture” and ‘‘Combined H:20 and organic matter,” 

4 Analyses also show 3.9 percent of CO; and 0.5 percent of Ci. 


Because large amounts of Ca(OH), are present in concrete, it 
might seem that the compositions given in table 4 should approximate 
those of the gels which coexist with Ca(OH)., as determined in this 
study. An examination of figure 7, however, shows that the re- 
ported compositions of the gels in concrete are notably and consist- 
ently lower in CaO and either higher or lower in Na,O than those 
along the boundary. 

A theory advanced recently by Hansen [10] to account for expan- 
sion of concretes containing “reactive” aggregates and alkalies also 
appears to explain why such low-lime alkali silicates may be present 
in concrete. This theory postulates that a membrane, permeable to 
sodium and hydroxy] ions and water but impermeable to silicate ions, 
is formed around the ‘‘reactive” aggregate particles. As a result of 
the transfer of the alkali hydroxide solution, an osmotic pressure 
great enough to cause expansion is produced. Also there occurs a 
reaction between the alkali hydroxide and the siliceous aggregate 
forming alkali silicates. 


VI. SUMMARY 


The only solids found, in the foregoing study of portions of the 
system soda-lime-silica water at 25° C, were Ca(OH), and a four- 
component gelatinous product of variable composition. In contact 
with crystalline Ca(OH), and solutions of increasing Na,O contents, 
this gel showed a composition varying between 0.003 Na,O:2.0 
CaO:1.0Si0,:2H,O (at 0.2 g of Na,O per liter) and 0.25Na,0:1.0 
CaO:1.0Si0O,:2H,O (at 20 to 101 g of Na,O per liter). By increasing 
the concentrations of SiO, in solution, at fixed Na,O levels (6, 13, 50, 
and 96 g per liter), the Na20:SiO, molar ratio of the gel was shown 
to increase to a value of about 0.2 and remain approximately constant 
~ — The CaO:SiO, molar ratio, however, decreased at all 

VacO levels. 
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The compositions of gels from concrete containing ‘‘reactive” 
aggregates are compared with those of the synthetic gels. 

It is shown that the low-lime gels formed in concrete approximate 
in composition those in the study in contact with solutions the Si0,: 
Na,O ratio of which approaches a value of 3. 
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LAMINAR FLOW AT THE INTERFACE OF TWO LIQUIDS 
By Garbis H. Keulegan 





ABSTRACT 


The velocity distribution in the laminar boundary layers at the interface of 
two liquids in relative rectilinear motion, the thickness of the layers, and the 
stress at the interface are determined. Numerical results are given for nine cases 
of liquids in contact, including identical liquids and liquids with varying degrees 
of dissimilarity in characteristics. The evaluation of the desired quantities is 
based on Prandtl’s boundary-layer theory, and is carried out by a method of suc- 
cessive approximations. The numerical results are those given by the second 
approximation. 





CONTENTS 
Page 
LBS ei sc eek ke ekg ce sel ck ce ass 303 
Td. erm Gante Gr GO POR a hn os enc oe ee cece 304 
III. Prandtl’s boundary-layer equation. -__...........--.-------------- 305 
IV. Approximate expressions for the stream functions. _.....-...---.--- 308 
Fe el a nec anmainaneeennwn 314 
VI. First approximate solution for the laminar boundary layers- - ----.-_- 317 
VII. Second approximate solution for the laminar boundary layers-_-_------ 318 
VIII. Quality of the approximations.....................-............. 319 
IX. Physically significant quantities of the layers._...........-.---.--- 321 
1. Local stress at the surface of separation.............--------.- 321 
a LO NY ee eRe 322 
3. Distribution of the velocities in the layers.............------- 325 
A. Leenevereal voiociien of tie Imyere_.... 3 -. i ec cl 326 
ks I Acti GALL Aiken h oeich ocbGhes oebed Sate i dercncnnonennse 327 


I. INTRODUCTION 


When a 4 ers of given density is flowing over a liquid of greater ™ 


density, the latter being at rest, three successive forms of the inter- 
face are discernible as the velocity of the moving liquid is increased. 
To take a concrete example, let us assume that the flow takes place 
in a closed rectangular channel with a lower pool; the upper liquid 
is fresh water and the lower liquid is a solution of salt or sugar in 
water. First, at small velocities, the surface of separation is sharp 
and distinct, indicating, since the indices of refraction of the liquids 
are different, that the densities are discontinuous at the interface 
and the flow at the interface is laminar. The flow in the central 
part of the upper liquid, on the other hand, may be laminar or turbu- 
lent, depending on the size of the channel. Secondly, if the velocity 
of the upper liquid is gradually increased, at some critical velocity 
the smoothness of the separating surface first disappears, and then the 
surface becomes covered with surface waves traveling in the direction 
of the current and with a velocity slightly smaller than the velocity 
of the upper current. When the waves first appear, they are long- 
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crested, that is, the length of the crests normal to the direction of 
motion is greater than the spacing between the crests. The crest 
lines are parabolic, concave upstream, and highest at the center of 
the channel, but flattening out toward the channel wall. At this 
stage of the flow—and this perhaps is significant—the waves are 
stable and show no tendency to break. The waves travel with 
practically no deformation, except the tendency to grow slightly 
larger during their motion. Finally, when the velocity of the upper 
current exceeds the critical velocity, the waves become sharp-crested 
and the crests shorter. The waves no longer are stable; that is, 
portions of the crests break away from the waves and are thrown into 
the upper current. These portions move forward and upward. In 
fact, this is the manner by which the mixing of the lower liquid with 
the moving liquid is brought about. When the velocity of the 
current is increased, still further within certain limits, the rate of 
mixing increases. At this stage of the flow, and this also is perhaps 
significant, the wave length is not affected by changes of velocity. 
These are the main results of observations made in a laboratory channel 
‘here the depth of the upper current is not over 10 centimeters and 
the length of the channel is not over 10 meters. 
A detailed explanation of these phenomena is lacking. A theoretical 
_/ approach supplementing experimental knowledge should prove help- 
: ful. A theoretical investigation should begin with a study of the 
stability of the interface on the basis of the actual conditions in the 
neighborhood of the interface. One of these conditions is a discon- 
tinuity of density at the surface of separation. Another is a pair of 
laminar boundary layers on the two sides of the surface of separation. 
In this paper the thicknesses of these boundary layers and the velocity 
distributions within them are determined by means of successive 
approximations, on the assumption that the interface is a smooth, 
plane surface. 


II. FORMULATION OF THE PROBLEM 


A liquid of indefinite height and having initially a uniform velocity, 
U, meets and flows over a still liquid of indefinite depth. In general, 
the two liquids have different physical characteristics; that 1s, their 
viscosities and densities are not the same. It will be assumed that 
the lower liquid has the greater density. If the pressure is everywhere 
hydrostatic, the interface or surface of separation is a horizontal 
plane. When U is sufficiently small the narrow region on both sides 
of the interface is a region of viscous laminar flow. Thus, the acti- 
vating forces at the point where the moving liquid comes first in con- 
tact with the lower liquid are tangential and act in the horizontal 
direction; there is no normal force acting on the interface and hence 
the interface remains horizontal. Every particle of liquid on the 
interface moves with a constant velocity, wu. The thickness of the 
viscous boundary layers in the two liquids is initially zero and increases 
with the length downstream (see fig. 1). In the upper layer the 
velocity in a vertical section increases from tu, to U; in the lower 
layer it decreases from %to0. Itis proposed to evaluate the variation 
of the thickness of the laminar layer with the length of the interface, 
the variation of velocity in the two layers, the shearing stress at the 
interface and the dependence of these quantities on the physical 
characteristics of the two liquids. 
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usU 
Laminar Boundary Layer U>uru, 
uzU 9 usu, Interface a 
xox 
Laminar Boundary Layer y! “u>u'>o 
'$ — 
v'<0,u'=0 

’ ‘ ° . oT enti 

Vea 'p lower liquid 


y’ V20,uU=0, y= 





| 
Figure 1 —Laminar boundary layers at the interface of two liquids. 


III. PRANDTL’S BOUNDARY-LAYER EQUATION 


In order to avoid negative values of ordinates two distinct sets of 
axes will be used, z, y for the upper liquid and 2’, y’ for the lower 
liquid. See figure 1. The z-axes of the two sets coincide and lie 
in the plane of separation of the two liquids, the positive sense being 
in the direction of U. The point at which the two liquids first come 
into contact will be taken as the origin. In the upper liquid, positive 
y is directed upward; in the lower liquid, positive y’ is directed down- 
ward. The density, the viscosity, and the kinematic viscosity of the 
upper liquid will be denoted by p, u and », respectively, and the 
corresponding properties of the lower liquid by p’, uv’, and v’, respec- 
tively. The parameter, r, defined as 


r=~V/[n'p’/upl, (1) 


is of special significance in this analysis, as will be seen subsequently. 

The basis of these computations is the boundary-layer equation 
of Prandtl for the two-dimensional flow of an incompressible fluid [1].' 
This equation with the law of hydrostatic pressure and the equation 
of continuity, when expressed in terms of the upper liquid, becomes 


ou, Ow Ou 1 Op 


Use Oy "dy p Oe’ (2 
x Soe . 
p dy + I= (3) 
Ou , Ow 
Bet dy” (4) 


where u and v are the components of velocity in the z and y direc- 
tions, respectively, and p is the hydrostatic pressure. 


1 The figures in brackets indicate the literature references at the end of this paper. 
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In view of the two-dimensional equation of continuity of an incom- 
pressible fluid, eq 4, the stream function, ¥, may be introduced: 


) re) 
=-3%, v= 5. (5) 


Since, in the present case Op/Oz=0, eq 2, in terms of the stream 
function, becomes 


Oy vy 40% ry dy 
Oy dzdy* Oa Oy" dy* 6) 
This can be changed into an ordinary differential equation by choos- 
ing appropriate new variables. The physical basis for the ordinary 
differential equation is the similarity of flow in two vertical sections. 
A necessary condition for the similarity of flow is that the velocity 
at the interface wu and the velocity at large values of y be constant. 
These conditions are fulfilled if the pressures are hydrostatic, which 
has already been assumed. Adopting a characteristic length 4, 
which is to be regarded as a function of x only, the new dimensionless 
variables for y and the stream function are 


n=ny/6, (7) 





and 


H(n)=—¥/U6, (8) 


where 7 is a dimensionless numerical constant as yet not specified. 
These variables are dimensionless. In terms of these, eq 6 now 
becomes 

ds ,,@H v@H 


Thus, if we select 6 to satisfy the relation 


5 dv 
a” T”’ 
or 
&=2n a (10) 
then, from eq 9, we have 
fH ?’H 
a3 +H ap (11) 


which is Blasius’ equation for the laminar boundary layer when the 
pressure is independent of z [2]. 

In eq 10 the dependence of 6 on 2 is given except for the numerical 
value of n. We determine this by considering the velocity at large 
values of y. Now, 


v= i aH 


oy dy’ 
hence if we specify that dH/dy=s1 for large values of 7, and this is a 


boundary condition for the upper fluid, we must then take n=1. 
These values we now adopt. 





_ 
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Summarizing, we have for the ee liquid 


L() =a tH GaH0 (12) 


where 


H=—fy n=4, (13) 


1% aw 


A similar set of equations may be established for the lower liquid by 
using primes to distinguish the hag 


and 


Li) =Ger +H! Gee =0) (15) 
where 
and 


ay)? to (17) 


The actual velocities and their first derivatives in terms of the 
dimensionless variables are 


ke 
=U GS (18) 
ou UH 
for the upper liquid; and 
fn ee i (20) 
Ow’ _ > on 
 ieual oy’ n’ “7.7 3? (21) 


for the lower liquid. 

For the simultaneous solution of the two basic differential equa- 
tions, eq 12 and 15, six boundary conditions are needed. These we 
now consider. In the upper liquid, wu equals U for y=~; in the lower 
liquid, u’ vanishes for y’=«. Accordingly, 


dH 
= a tig (22) 


and 


ee ~ ere: (23) 
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The stream function, ¥, may be interpreted as the flux per unit width 
between a point y of the upper liquid and the interface; and y’ as the 
flux per unit width between a point y’ of the lower liquid and the 
interface. As we have supposed that the interface is the horizontal 
plane y=y’ =0, we must have 


H=0, n=0; (24) 
and 
H’ = 0, 9=0. (25) 
The velocities are continuous at the interface; hence 
oi dH’ 
~ dy’ 37 ’=(), (26) 


Again, the shearing nila are continuous at the interface; hence 
Ou ou’ j 
ewe 
u @H u’ @H’ 


, 
—— Gwe ees eee SO == =f) 
Vv dn? Vv’ dn’ 2 27=7 ’ 
or, using eq 1, 
eH. @H _, ae 
de "dy® n= =(), (27) 


These are the six boundary conditions of the problem. 


IV. APPROXIMATE EXPRESSIONS FOR THE STREAM 
FUNCTIONS 


To solve the differential equations, eq 12 and 15, we shall resort to 
a method of approximations which is a modification of the method 
used by Pohlhausen for the solution of the Blasius plate problem [3]. 

In this method certain approximate boundary conditions for large 
values of 7 and 7’ are required. First, consider the upper moving 
liquid. Select the value n, of 7 such that 


] -7 s ~ €, if n= Ns, (28) 


where e is a small quantity, say 1/100 or any other smaller fraction. 
Denote the value of H corresponding to », by H,. Now, writing 
eq 12 in the form 
fH 1dH 
dy H di’ 


and integrating between n=7, and n=, we have 


-(¢ a) =— ae GH 
* , H die 
Integrating the right-hand member Eee parts, and neglecting the 


term multiplied e, and making an obvious substitution from eq 29, 
using eq 28, we obtain 


d 1 (@ -1 &H 
(Ge) Hae), + J, He dy 7 


(29) 


Co 
suc 


Ob 


} 


lar 
def 
beg 
cou 


fro) 
ysis 





os 


J 


Flow at Liquid Interfaces 309 


Treating the integral in the right-hand member again in the same 
manner, we obtain 


dH\ _(fH\(1_ 1 21 PH 
1-(&), =(Ge) (az a)-3{. He dy dn. 


If ¢ is sufficiently small, that is, if H, is sufficiently large, the integral 
in the latter equation may be discarded, leaving 


(8) NG) 
dn Jy \da j»AH, He (30) 


This is the approximate relation which exists between the first and 
second derivatives of H when n= 7. 

If we now suppose that H, is so large that 1 can be neglected com- 
pared to the square of H,, then, in view of eq 28, 29, and 30, we must 
have 





dH : 

“dn =1—e, n=MNsy 

?H 

dg =e N=MNsy? (31) 
PH . * 

dy? =—6H,’, i 


Consider now the lower still liquid. We select the value n,’ of 7’, 
such that 


* hp aes 
‘dr’ Se’, if n'>n,’, (32) 


where e’ has the same general meaning as «. Denote the value of 
H' corresponding to 7,’ by H,’. For values of n’ larger/than n,’, eq 
15 may be written in the form 


PH’ __ yy PH’, 
dy’* ae 8 dn’? 
Obviously, then 
dH’ ? ) 
i” * Bate gr Fes 
d’H’ , , / , 
Gq? eh, =} (33) 
ere Ee” 
dnt =e", n =M4s - 





Now, in the method employed here, the end conditions for the 
large values of 7 and 7’ will be taken from eq 31 and 33 after assigning 
definite values to « and e’. The values assigned to « and ¢’ have a 
bearing on the accuracy of the approximation and also govern the 
course of the computations. If « and ¢’ are assigned values different 
from zero, the accuracy of the approximation is increased, but the anal- 
ysis becomes more complicated. In order to simplify computations, we 
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shall take e=e’=0 and adopt the following end conditions. For the 
upper liquid, we have 
dH ix @?’H 
dn ’ dn =0, ae 0, n=, (34) 
and for the lower liquid, 
dH’ d°H’ dH’ Pe 4 
ar dn? =0, or »1 = - (35) 
These now replace the conditions implied by eq 22 and 23. With 
this choice of the end conditions, y, and 7,’ have simple geometric 


interpretations. Let y, and y,’ be the values of the distance from the 
interface corresponding to 7, and »,’, respectively; that is, let 


Ys=Ns V (2v2/U), (36) 
and 
ys’ =n, ¥(2v’2/U), (37) 
then y, and y,’ are the thicknesses of the two laminar boundary 
layers at the interface. 


We shall also need the following end conditions from eq 24, 25, 26, 
for the upper liquid, 


dH 


H=0, a ™ n=0, (38) 
and for the lower liquid, 
H’=0, a =a)’, n'=0, (39) 


where 
p= Ay’ =Up/ U. 


We are now ready to develop the general expressions for the stream 
functions H and H’. Since from eq 15 and 34, d*H/dy* vanishes at 

=0 and n=7,, then it must be possible to express this function of H 
as a summation of polynomials where each polynomial vanishes for 
n=0 and n=7;. Suitable polynomials may be selected in various 
ways. In order to extend Pohlhausen’s computations [3], we select 
the following oscillating algebraic functions as the expansion poly-| 
nomials for this analysis: 


N,=— a(a—1)=a—a’, 
N2=+a(a—1)(2a—1)=2a°—3e’+a, 
N;=—a(a—1) (83a—1)(3a—2) =2a—11a*+ 18a--9a', 
and, in general, 
N,=(—1)"a(a—1) (na—1)(na—2)(na—3)( . . . )(na—n+1), 


where 
a=n/n-. 








(35) 
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Using only the first three functions, we put 
a “(—6N,+32k,N;—30kN,). (40) 


Here a, 7s, k2, and k; are constants to be determined. Integrating 
once With respect to 7, we have 


Fatt kletkwh) (41) 


where J,, J. and J; are functions of a: 

J,;=1—307+ 20%, 

J,= 16a°— 32a°+ 16a4, 

J3;=1—3007+ 110a®—135a*+ 540°. 
It must be remarked ‘that J; is derived from N, by integrating with 
respect to » and selecting the constant of integration so that J;,=0 
for a=1, and similarly for J; and J3. Equation 41 thus ay igen 


satisfies the second end conditions in eq 34. Integrating eq 41 wit 
respect to n, we have 


a= dot ain Gr + keGrt kaGh), (42) 


where G,, Gj, and G; are functions of a: 


Gi;=a— a? +50", 


y= a 8a + of, 


G3=a— 10a?-+ >? a! —27a°+ 9a. 


In this case, G, is derived from J, by integrating the latter with re- 
spect to 7 and selecting the constant of integration so that G,;=0 
when 7=0, and similarly for G, and G3. The constant Gy is intro- 
duced in order to satisfy the second end conditions in eq 38. Inte- 
grating eq 42 with respect to 7, we have 


H=a9nsFo+ Qn (P+ keF + keFs), (43) 
where Po, F;, F2, and F; are functions of a: 
F\=a, 
1 1 1 
Py=5e—Ge+79% 


. 8 8 
2 =g0'— alt is” 
5 


Fy=5o? 


cere sec Ta 
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where Ff is derived from G, by integrating the latter with respect to 
n and selecting the constant of integration so that F,=0 when »=0. 
and similarly for Fy, F,, and F;. It is obvious that the first of the 
end conditions in eq 38 is satisfied. The functions N, J, G, and F 
are tabulated in tables 1, 2, 3, and 4. 


TABLE 1.—Values of the N functions 














ioe 4 oe | Nila) X10 | Nala) X10 
0.00 0. 0000 0. 0000 0. 0000 
= | 0475 14275 , 7030 
Lo 0900 : 7200 1.0710 
15 . 1275 . 0925 1. 0869 
fear. | : 1600 : 0. 8960 
| 125 ‘1875 9375 ‘5859 
:30 "2100 "8400 2310 
eee | 2275 6825 | —.1081 
40 2400 -4800 | —.3840 
145 "2475 2475 —. 5631 
‘50 ‘2500 ~, | —.6250 
55 12475 —.2475 | | —. 5631 
60 | 2400 —! 4800 —. 3340 
65 2975 —.6825 | —.1081 
‘70 : 2100 —.8400 | 2310 
75 | 1875 —,0875 | 5859 
| 80 -1600 | = —. 9600 8960 
\. So 1275 —, 8925 1. 0869 
| 90 0900 | —.7200 | 10710 
ae | 10475 —.4275 | 0.7030 
| 1,00 .0000 | —.0000 0000 











TaBLeE 2.— Values of the J functions 




















i | | 
a Ji(a) | J(a) J3(a) 
0.00 1.0000 | 0.0000 1.0000 | 
05 0. 9929 . 0344 0.9451 | 
10 . 9720 | . 1296 . 7971 
15 . 9393 . 2632 6312 | 
| 2 «| 8960 | . 4096 4842 | 
| m4 . 8437 . 5624 3600 | 
30 7840 | . 7056 3075 | 
35 .7183 | . 8280 .2979 | 
40 . 9216 . 3366 
| 45 5747 . 9800 . 4105 
50 1. 0000 . 5001 
55 4253 | 0.9800 . 5904 
60 3520 | . 9216 . 6633 
65 2317s . 8280 . 7008 
70 +} 2160 . 7056 . 6942 
| 7 1563 . 5624 . 6309 
80 . 1040 - 4096 . 5280 
| 85 "0607 2032 | 13675 
90 . 0280 . 1296 . 2037 
9 | . 0073 0344 | . 0615 
| 1.00 . 0000 - 0000 . 0000 











*t to 
=0, 
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TaBLE 3.—Values of the G functions 











a | Ge Gs (a) Gia) | 
0.00 0. 0000 0.0000 0.0000 | 
05 10499 ‘0011 10480 
"10 0990 “0045 ‘0927 | 
15 ‘1469 0144 1270 
20 ‘1928 "0299 "1569 
- 25 . 2363 . 0552 . 1760 | 
- 30 - 2770 . 0869 . 1942 j 
35 3146 "1254 ‘2093 
40 3488 . 1691 2255 
45 3793 . 2168 2444 
50 4062 "2669 2641 
55 4204 "3165 2932 
60 4488 ‘3640 3244 
65 4077 3500 | 
70 4771 Y .4464 3937 
75 4863 ‘4781 4272 
80 4928 ‘5024 4559 
85 4969 . 5191 4783 
90 4990 * 5288 4925 
95 4999 ‘5326 4999 
1, 00 5000 "5333 5000 




















TABLE 4.—Values of the F functions 











a F; (a) F2 (a) Fs (a) 
0.00 0. 0000 0. 0000 0.0000 | 
05 0012 . 0000 aa 
.10 . 0050 . 0001 .0047 | 
15 0112 . 0006 0106s 
20 0196 . 0016 0176 
25 0303 . 0037 0261 :~«—s| 
30 0432 . 0073 0351 
35 .0127 0449 
40 0746 . 0200 0556 | 
45 0928 . 0297 0692. | 
50 1125 0415 0807 
55 1334 . 0562 0995 
60 1554 . 0733 1096 
65 1782 . 0925 1267 
70 2018 . 1138 1457 
75 2258 .1370 ° 1661 
80 2505 . 1618 1881 
85 2753 . 1874 2118 
90 3001 . 2134 2357 
95 3750 . 2398 2608 
1, 00 3500 . 2666 2857 




















In these expressions there appear two types of constants. The 
constants k, and k, are of one type and may be referred to as the ex- 
pansion constants, since these are introduced because of Nz and Ns. 
The constants dp, a, and 7, are of the other type and may be referred 
to as the basic constants, since dp relates to the interfacial velocity 
U, A, to the local shear and 7, to the thickness of the laminar boundary 
layer in the upper liquid. 

Reasoning in a similar manner, we have for the lower, still liquid, 


oon 4, (—6N,+32k2!N,—30k,’N;), (44) 
Gey thy Tethys), (45) 
dH’ . 
oS =’ + y's’ (Gi +key’ Ga+k,’ 13) 5 (46) 


and 
H’ =p’ ny’ Foy’ 5"? Fi + hee’ Fo +h’ Fs), (47) 
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where N, J, G, and F are the same functions as before but now depend 
on the independent variable of a’=7’/n,’. Again, kp’ and k;’ are the 
expansion constants, and a)’, a;’, and 7,’ are the basic constants for 
the laminar boundary layer in the lower, still liquid. 

Equations 43 and 47 are the approximate expressions for the stream 
functions H and H’. The unknown elements are the basic and expan- 
sion constants, which we now proceed to determine. 


V. EVALUATION OF THE CONSTANTS 


There are a total of six basic constants and four expansion constants 
in the expressions for H and H’, and thus 10 relations are required for 
their determination. Four such relations are obtained immediately 
by considering the first of the end conditions in eq 34 and 35 and the 
continuity conditions, eq 26 and 27. These give, first, 


Ao+ainsAyp= 1, (48) 
where 
16 ; 
second, 
Qo’ +,'n; Ay’ =0, (50) 
where 
Ad =3(1+35 be’ +h’) (51) 
third, 
A=) ; (52) 
and fourth, 
a,=—ra,’M; (53) 
where 
M=(1+k,’)/(1+ks),  1?=(u'p’)/(up). (54) 


The remaining six relations may be obtained from eq 12 and 15 by 
using one of two methods. 

In the first method the left-hand members of eq 12 and 15 are mul- 
tiplied successively by Ni, No, and N;, and integrated with respect to 
a or a’ between the limite 0 and 1 and set equal to zero, since L is 
identically equal to zero. That is, putting 


PH. ?H : 
L= dp +2 dr’ (55) 
and 
_@H’ @H’ 


(56) 
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the desired six relations are 


S3LN,da=0, (i=1, 2, 3); (57, 58, 59) 
and 
1 L’N, da’ =0, (i=1, 2, 3). (60, 61, 62) 


A simpler method would be to use Q;, Qo, and Qs instead of Ni, No, 
and N3, where 


Q=+1, OSaora’Sl1; 
Q.=+1, OSaor a’S1/2; 
Q.=—1, 1/2Saor a’s1; 
Q,=+1, OSaor a’S1/3; 


Q,=—1, 1/3Saor a’ $2/3; 
Q@=+1, 2/3Saor a’S1. 


The second method, which is the more familiar one, utilizes the 
following six relations [4]: 


JS Ida=0; (63) 
5 L’da=0; (64) 
dL 
a n=Ms} (65) 
dL 
dL’ 
ay" n=" ; (67) 
dL’ 
G0 1'=0. (68) 


The constants of the present problem have been determined, by 
using both of the above methods, and the results agreed to the same 
order of accuracy. The résults obtained from eg 63 to 68 will be 
given here. 

From eq 40, 41, and 43, we find from eq 55, 


L(n) = (—6N,+32k:N)—30k,N5) 


+aoains (Si Fo+ koJoFo+ kyJ3Fo) 
+ ayn? (SF; + heoldJoF + Feds] +klJsF; + FJ;]) 
+a;"n;? (ke? I oF, + kks[J2F'3+ FJ3)+ k?J3F 3) . 


Substituting in eq 63 and integrating, we find that 


y?== Ajdo0;"n,’ + Apa,?n,’, (69) 
where 


A= (Ao. + Ao2k2+ Aogks)/(1+ks), (70) 
A= (Ay, + dak. + ayaks+ dook2? + degkoks + Asks”) /(1 +s). (71) 


and 
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The values of the coefficients with the double subscripts occurring 
in eq 70 and 71 are 


d=0.15000, a=0.029365,  a_=0.0308, 
2 =0.26667, dy2=0.07619, og = 0.0758, 
Og=0.21425, dy3= 0.07625, dz3= 0.0400. 


By treating eq 64 in a similar manner, we find that 


a," = Ay'do'ay"n,” + Ay’ a,”"9,”» (72) 
where 
Aj! = (dort doko! +dogks’)/(1+hs’), (73) 
and 


Ag! = (Ay: + ike’ + digks’ + ook!” + Angle! key’ +assks”")/(1 +k’). (74) 


Since from eq 55 


dL dH dH@H , ,,@H 
dy dat | dy dy + dy 


and 
N=, ON? 1, oN) =2, n=0, a=0, 
M1, one 1, oi —2, 1=", a=1, 
we find from eq 65 
= (6+32k:+ 60k) =0, (75) 
and from eq 66 
“1(—6-+32k,— 60ks) = —dga (1+ hy). (76) 
Adding eq 75 to eq 76, we have 
64k. = —don,?(1+ks). (77) 


Since », is finite and a; does not vanish, we have, from eq 75, 


We have from eq 77, using eq 48, the relation 


eee 2 
64k, = — fae] +h) /(1 + 48k +h) } (79) 


where 7, does not appear. It may be emphasized that eq 78 and 79 
are two relations expressing the expansion constants k, and ks in 
terms of the basic constants a and 4. 
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In a similar manner we obtain from eq 67 and 68 the two relations 


60k,’ = —6—32k,’, (80) 
4a,’3 2 
oak,’ = — Se] (a +hy’)(1+ Fp’ +h’) | (81) 


Equations 48, 50, 52, 53, 69, 72, 78, 79, 80, and 81 are the desired 10 
relations for determining the 10 constants appearing in the expressions 
for H and H’ given by eq 43 and 47, respectively. In principle the 
determination of the constants in a direct manner is possible, but 
since this is tedious, it may be carried out by successive approxima- 
tions. The successive approximations to the values of the constants 
give the successive approximations to the present boundary-layer 
problem. 


VI. FIRST APPROXIMATE SOLUTION FOR THE LAMINAR 
BOUNDARY LAYERS 


In the first approximate solution it is assumed that the expansion 
constants kh», ks, ko’, and k,’ are zero. 
Consider the relations 


Ay+4\n,Ao=1, (48a) 

Ap’ + y'n,’ Ao’ =0, (50a) 

Qyp= Ay’, (52a) 

a,=—ra,'M, [r’=(u'p’)/(up)], (53a) 

0? = AjQo0;"n,’ + A2a,*n,’, (69a) 

%: ay’? = Ay’ do'd,'?n,/? + Ag’ ay*n,”*. (72a) 


Here the constants Ao, Ao’, M, A;, A;’, A; and A,’ are functions of kp, 
ls, ko’, and k;’ and may be determined from eq 49, 51, 54, 70, 71, 73, 
and 74. For ke, ks, ke’, and ks’ all equal to zero, we obtain 


Ap= Ap’ =0.50000, 
A,=A,’ =0.15000, 
A,= A,’ =0.023936, 
M=1. 


Now, the determination of the unknown quantities a, a, 7, and 
ty’, dy’, n»’ from the above six relations, 47a to 72a, is an easy task 
once a specific value is assigned to r. Separate determinations were 
made by selecting nine values for r, namely, r?=0.00, 0.01, 0.10, 0.3162, 
1, 3.162, 10.0, 100, and oo. The results are shown in table 5. 
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TaBLE 5.—First approximate values of the constants ao, a), a’, 9, and n,’ 
ri=(y'p’)/(up) 





























ri a1 a’, ao or ao’ | Ne mn! | A : ed 

| ieee 0.6038 | 1.000 |_....._... ee eee 0.995 | 
01 0.0540 0.9285; 261i | 3.433| 2552 959 
10 1388 4300| .8085| 2.744] 3.678| 2.582 893 
3162 2031 3614| :7100| 2845] 3.923) 2.505 838 
10 2727 2724| :5910| 2995] 4.331 2. 616 772 
3. 162 3397 1911 4646 | 3.146; 4851| 2630 679 
1 3916 1238} .3479| 3.824| 6.605| 2.672 590 
100 4494 0449 :1773| 3.655] 7.844, 2.750 423 
00 4ss0 | ..._.- "000 06 fF ci... | 2.884| _._... 

i 








The expressions for H and H’ from eq 43 and 47, respectively, 
wherein the constants k and k’ are put equal to zero and the values 
of the basic constants are taken from table 5, constitute the first 
approximate solution of the laminar boundary layers at the interface 
of two liquids, one of which is still and the other moving. More 
specifically, in the upper liquid 


dH 
H= agnsl" (n/ns) +a,'n.F; (n/ns); dy tan (n/ns), 
and in the lower liquid 


, 


, , , , / 72 , ul , / dH P4 / , , / 
Hl’ =o" ns’ Fo(n' /ns’) +-0n/?ns"?F s(n’ /ns'), Gy” 0 tat G;(n’/n,’). 


The numerical values of Fy, F;, and G; are given in tables 3 and 4. 

It may be of interest to note that the case r= can be physically 
interpreted as representing the case in which the lower liquid has 
an infinite viscosity, which is obviously the case for a liquid flowing 
over a plate. The solution obtained here for r= © is identical with 
Pohlhausen’s fourth approximation. In the present analysis this con- 
stitutes only the first approximate solution. 


VII. SECOND APPROXIMATE SOLUTION OF THE LAMINAR 
BOUNDARY LAYERS 


Inserting the first approximate values of do, a, and do’, a,’ from 
table 5 in eq 78, 79, 80, and 81, we obtain a set of values for kz, h;, 
k,’, and k,’, after discarding the squares and the higher powers of 
these constants. The results are given in table 6 for the various 
values of r?. 


TaBLe 6.—Second approximate values of the expansion constants ks, ks, kp’, and hi’ 

















rts (yup) /(up) 
| | | | 

Am 4 Oe chk te) bight ky «| OM 
ot Baer a Lier See ae —0.2600 | +0.0382 |............ 
101 | —0.1354 | 010278 | 22800 | +.0382 |” 1.0681 
10 | —.1259 | 0328 | 2600 | +.0382 | 1.0787 
13162 | —.1162 | —.0380 | —.2600 | +.0382 | 1.0816 
10 | —.1037 | —.0446 | —.2800 | +.0382 | 1.0869 
3.162 | —.0807 | —.0522 | ~.2600 | +.0382 | 1.0956 
10 | =.0734 | —.0608 | —.2600 | +.0382 | 1.1057 
100 | —0429 | -.0771 | —.2600 | +.0382 | 1.1252 
eo | 0000 | = 2000 | nn. | a... F cdoierts 
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With these, corresponding values of the coefficients A), A’, Aj, 
A;’, Ae, Ao’, and M are computed. The results are shown in table 7. 
With the values of the coefficients thus known, the second approxi- 
mate values of the basic constants are determined by the same 
equations as for the first approximate solution, that is, eq 48a to 72a, 
inclusive. The results of this determination are given in table 8. 

Expressions for H and H’ from eq 43 and 47, respectively, where 
the expansion constants are given in table 6 and the basic constants 
are given in table 8, constitute the second approximate solution of 
the laminar layers boundary problem for the interface between two 
liquids. 


TABLE 7.—Second approximate values of the auailiary constants Ao, A,, Az, Ag’, Ax’, 

















ely, and A,’ 
tues r= (y'p!) (up) 
irst 
face | 4 | , ’ , 
lore r Ap | Aj A: Ag A, A: | 
| | 
fe ers eer pee Memb pe ne 0.3808 | 0.08562 | 0.01336 | 
ee 0.4139 | 0.1110 | 0.01830 . 3806 ‘08562 01336 | 
ae 4164 1131 ‘01871 3806 ‘08562 .01336 | 
| 3162 4190 1155 .01917 . 3806 08562 01336 | 
| 1.00 4224 1181 01970 . 3806 . 08562 :01336 | 
| 3. 162 . 4260 , 1212 . 02038 3806 08562 01336 | 
10 4304 1250 .02107 | .3806 08562 01336 
100 4382 1321 02245 3806 08562 01336 
, 00 4600 1427 yo Ppdaaes cs Bota Ee meant 
: | 


rom 


) ks, 


3 ol 
ious 


1 ky! 



























TaBLE 8.—Second approximate values of the basic constants ao, a), a;’, ns, and no 


























r2= (u'p’)/(up) 
r aq a,’ Qo Or ao’ Ne n A, | H, 
| 
| 0.00 |...-..---- —0.5906 | 1.0000 |_..-..-.-- 4. 449 3.012 1.041 | 
01 0. 0567 —.5329 | 0.9282 3. 069 4.576 3.004 0. 993 
10 . 1456 —. 4287 , 8082 3. 163 4.953 2.991 1936 | 
3162 . 2144 —. 3527 7088 3, 241 5, 280 2. 990 875 | 
1.0 . 2899 ~ ; 3. 358 5.793 2. 991 . 798 
3. 162 . 3611 —. 1854 : 4619 3. 498 6. 545 2. 996 ; 707 
10 . 4174 —.1194 3446 3. 648 7. 582 3.016 . 602 
100 , 4817 —. 0428 .1739 3.913 | 10.672 3. 055 434 
00 te Reps Se . 0000 “cS eel: Sani: seen 














VIII. QUALITY OF THE APPROXIMATIONS 


The quality of the approximations may be studied conveniently by 
substituting the approximate expressions for H and H’ and for their 
derivatives into the differential equations, eq 12 and 15, and then 
considering the magnitude of the remainders. Let the remainders 
be denoted by AL and AL’: 


_@H 2H 
AL ~~ dn +H 
and 
es 2 
AL =q,0 tH dn” ° 
If the ratios 
AL/(@H/dny?) max and AL’ /(d?H’ /dn’*) wax 
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become smaller in the successive approximations, the magnitude of 
decrease is a measure of the improvements in the approximations, 
We have computed these remainders only for the two extreme 


0.24 a oe 
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‘in 3 ae \ 


1St app. \Y | 
4L | 
Behe ee, 


-0.06 a tak: | a sd 4 
ie) 


| 2 3 4 5 


Figure 2—Remainders of the first and of the second approximations for the case 
r#=00, 
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Ficure 3—Remainders of the first and of the second approximations for the case 
ré=0, 
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cases represented by r?= and r’=0. The values of the remainders 
together with the values of the third derivatives of H and H’, are 
given in figures 2 and 3, respectively. The case of a liquid of finite 
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viscosity moving over a liquid of infinite viscosity is represented b 
r=. As mentioned before, this is mathematically identical wit 
the Blasius problem of flow over a plate. In examining the curves 
in figure 2, we see that the yo. oot rer of the oscillation of the remainder 
in the second approximate solution is only about half as large as in the 
first approximation, indicating that the improvement in the second 
approximate solution is appreciable. In the second approximate 
solution, the maximum value of the remainder is 0.027; the maximum 
value of the third derivative, 0.25. Itis thus seen that the remainders 
are relatively small. In this same figure we give also the curve of the 
third derivative of H from the exact solution by Dryden [5]. The 
agreement between the second approximate solution and the exact 
solution is rather close. The velocity gradient at y=0, from the 
second approximate solution, is 


1 ou 

U oy 
In the exact solution the corresponding coefficient has the value 
0.3320, thus indicating that the error in the value of this quantity 
from the second approximate solution is only 0.39 percent. 

A liquid of infinite viscosity moving over a liquid of finite viscosity 
is represented by r?=0. Physically, it would seem that this is the 
case of a solid surface moving in a still liquid. The curves in figure 3 
also show that in this case the remainder of the second approxi- 
mate solution has been reduced and has become about half as large as 
the corresponding quantity in the first approximate solution. But the 
accuracy obtained is not as satisfactory as in the case r?= ©, since 
the remainders are relatively larger. Apparently to obtain the same 
degree of accuracy as was obtained in the case r?= ©, and with the 


second approximation it will be necessary to introduce the expansion 
functions k,N, and k;N;. 


IX. PHYSICALLY SIGNIFICANT QUANTITIES OF THE 
LAYERS 


=0.3307 VU/(vz). 


Among the quantities of the laminar layers at the interface of the 
two liquids, one at rest and the other moving over it, the following 
are of practical significance: (1) The intensity of the viscous stresses 
at the surface of separation, (2) the thickness of the layers, and (3) 
the distribution of the velocities in the layers. The computation of 
these quantities will be made, using the second approximate solution 
of the analysis. 


1. LOCAL STRESS AT THE SURFACE OF SEPARATION 


The stress in question is determined from the velocity gradient of 
the upper liquid at the surface of separation: 


o (Ge 
T=p 5 dr? 4 


From eq 19 
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or 
2»/d°H\ p Ty? 
T= TH Ts ha Us 
6U dn’ 02 


Substituting from eq 14, 
_ pf ©H\ (Uz\-*p 5» 
a . x( dn? ), \ v ) 2 U 


If we use the velocity gradient of the lower liquid, we obtain in the 


same manner, 
ea d’H’ Uz\~%p’ rp? 
sme a =, 5 tom A Bi 


Putting 
5 HH Ps 
\ (Fa Je V¥2a,(1+ks) =s; 
and 
PH’ - ’ 
.dn’" /o 


we have the local stresses 


2 
T=39 *) £ U?, ($2) 
v 2 


ol CRN TE eee ’ 
roe ( =) 5) U?. (83) 


The computed values of s and s’ as functions of r? are given in 
table 9. 





or 





TABLE 9.— Values of the constants s and s’ in the expressions of the local stress 
ai the interface, eq 82 and 83 


r= (n'o")/(up) 











r? 8 | 8 

0. 00 ie sill 0. 8672 

. 01 | 00782 | . 7825 

.10 .1992 | . 6295 

.3162 | : ee . 5178 
1.0 3017, .3918 | 
3. 162 .4840 | "Bee 
10 .6543 | 1788 =| 
100 .6287 | .0628 | 
oo) 6612 . 


2. THICKNESSES OF THE LAYERS 





Since the velocities in the laminar boundary layers approach the 
limiting values in an asymptotic manner, the definition of the thick- 
ness of the layers becomes a matter of convention. Of particular 
utility are the following definitions. See figure 4. 


2 j, ” (U—u)dy=(U—u)&i, (84) 
0 
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Fiaure 4.—The gradient thicknesses of the interfacial laminar boundary layers. 
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and 

v.’ 

2 au’ dy’ =Uby', (85) @ 

the 
the first applying to the upper moving liquid and the second to the § u/I 
lower still liquid. It will be appropriate to speak of 6, and of 6,’ as res 
the gradient thicknesses. Like the displacement thicknesses, the 
values of 5, and 6,’ would not be affected very much by the approxi- 97" 
ee and this obviously, can not be said of y, and y,’. Using the 
relations 





y=n, y’=1's', ‘i 
6=y2n]U, 8’ = y2v2]U, 
u=UdH/dn, u’ = UdH' /dn’, ; 
we find from eq 84 and 85 that | 
— p—H, va 
a=sys 1—a VU’ 
and 50a 
H,’ |/v'z 
ee, ee ee 
61 2y2 do U 
Putting 
=—9 n.—H, Pee 5H,’ T 
n=2 2 1—a& 2/2 do ’ Al 


we have for the gradient thicknesses, 
=n and 6,’=n’ ia (86) 


The values of n and n’ are given in table 10 as functions of r?. We 
obtain for the boundary thicknesses, introducing the expressions 


m=-+2n., m! = J2n,’: 
yam], y; =m’ ea (87) 


The values of the coefficients m and m’ as functions of r? are given in 
table 10. 


TABLE 10.— Values of the constants n and n’ in the expressions of the gradient thick- 
ness and m and m’ in the expressions of the boundary thicknesses of the two lam- 
tnar layers, eq 86 and 87 














r2=(y'p’)/(up) 
ri n | n’ m m’ 
| 
| 
BE ee: Bae ah Sale 2 6. 291 
01 2. 270 3.026 4, 339 6. 470 
10 2. 344 3. 304 4. 472 7.003 
3162 2. 428 3. 493 4, 582 7, 466 
1.0 2. 531 3.831 4.748 8. 191 ; 
3. 162 2. 664 4. 328 4.946 9. 255 Ine 
10 2.813 4.940 5.158 | 10.721 
100 3.071 7. 063 5.533 | 15.090 
© i es eae SOM 1 Zu. 
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3. DISTRIBUTION OF THE VELOCITIES IN THE LAYERS 


These velocities are determined from eq 42 and 46. The results of 
the computations for the upper liquid are given in table 11, where 
u/D is tabulated as a function of y/y, for various values of r?. Similar 
results for the lower liquid are given in table 12. 


TaBLE 11.—Velocity distribution in the interfacial laminar layer of the upper, 
moving liquid 


The boundary thickness y, is given by eq 87 


















































r2= (u'p’)/(up) 

ee ya 

lv. ) 0.010 0. 100 0. 316 1.000 3. 16 10 100 © 

u/U 

| ice , eT LIE De eos ore = 

| 00 0.9282 | 0.8082 | 0.7088 | 0.5888 | 0.4619 | 0.3446 | 0.1739 | 0.0000 

a . 9449 - 8521 7747 . 6807 . 5808 . 4862 . 3466 . 1994 

| 2 . 9602 . 8929 . 8363 . 7667 . 6917 . 6204 . 5122 . 3935 
3 . 9735 . 9279 . 8892 . 8412 . 7891 . 7387 . 6609 . 5723 
.4 . 9838 . 9557 . 9316 . 9015 . 8685 . 8359 . 7850 7271 
5 . 9913 9758 . 9626 . 9457 . 9273 . 9089 . 8797 . 8440 
.6 . 9961 . 9889 . 9826 . 9747 . 9660 . 9573 . 9432 . 9252 
i . 9988 . 9962 . 9933 . 9910 . 9878 . 9848 . 9801 . 9727 
8 . 9998 . 9991 . 9985 . 9980 . 9973 . 9966 . 9960 . 9938 
9 | , 9999 . 9999 . 9998 . 9998 . 9997 . 9996 . 9996 . 9995 
10 | 1.0000 | 1.0000 | 1.0000 | 1.0000 1. 0000 1.0000 | 1.0000 | 1.0000 














TaBLE 12.—Velocity distribution in the interfacial laminar layer of the lower, 
still liquid 


The boundary thickness y,’ is given by eq 87 












































r= (u’p’)/(up) 
| 
ly mr 0.00 | 0.010 | 0.100 | 0.316 | 1.000 | 3.16 10 100 
| | 
w/U 
0.0 1.0000 | 0.9282 ] 0.8082 | 0.7088 | 0.5888 | 0.4619 | 0.3446 | 0.1730 
cf 7337 | .6811 | .5030 | .5200 | .4320 | .33899 | .2527 | 1276 
3 4980 | .4624 | 14023 | .3530 | 2932 | ‘2301 1716 | .0866 
3 3121 2307 | 2524 | 12213 | 11838 | :1441 1075 | :0542 
4 1765 1639 | 11427 | .1250 | :1030 | 0814 | :0608 
5 ogss | .0822 | 0717 | :0628 | 10522 | ‘0404 | ‘0305 | ‘0154 
6 0369 | :0343 | :0208 | :0262 | ‘0217 | ‘o171 | ‘oust 
7 0119 | .0112 | :0098 | 0092 | ‘0070 | ‘0055 | ‘041 | ‘0020 
8 0026 | 0025 | .0022 | :o021 | ‘016 | ‘oo11 0004 
9 0007 | .0006 | .0005 | :0004 | :0003 | :0002 | ‘0001 | ‘0001 
1.0 0000 | :0000 | :0000 | 0000 | 0000 | ‘0000 | ‘0000 | ‘0000 





A more suggestive picture of the velocity distribution is obtained 
by plotting u/U against y/y,, y now being the distance from the surface 
of separation to a point, above or below the surface, and y, being the 
boundary thickness of the layer situated in the upper liquid. Figure 
) gives the velocities for identical liquids, uy’ and p=p’. Figure 6 
gives the velocities for v’=7.85v and p’=1.14p. This case will be 
realized if we take water for the upper moving liquid and aqueous 
solution of 56 percent glycerol for the lower, still liquid, both being at 
20°C. These two figures demonstrate the very large influence that the 
increased viscosity of the lower liquid has on the velocity distribution. 
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Ficure 5.—Interfacial laminar velocity distribution for the case r?= 1, 
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Ficure 6.—Interfacial laminar velocity distribution for the case r?= 10. 


X. TRANSVERSAL VELOCITIES OF THE LAYERS 


The expressions of the transversal velocities, v and v’, in terms o! 
H and n, may be derived conveniently from the condition of continuity 
eq 4. Since at the interface, v and v’ vanish, we have 


Y Ou 
o=— Sx UW: 








of 
ity 
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This yields, using eq 13, 14, and 18, the relation 
1 ?H 
ser =f dy ndn. 
Integrating by parts and noting that H=0 when n=0, we obtain 


iv dH 
a ei! (88) 


which is the expression for the transversal velocities at the point 
in, 6). 
The corresponding expression for the lower fluid is 


5’0’ v __ dH’ 4 

~~ dn ~yq — EH’. (89) 
Thus, at the points of the lower liquid and in the region outside of the 
layer, »’>»,’, there is a weak current moving normally upward 
toward the interface. The strength of the current is given by 


helmets col (90) 


In a sense, the boundary layer of the upper liquid acts as a pump, 
Swe the small portions of the lower liquid to the level of the inter- 

face, then causing these portions to move horizontally. So long as 
these motions as required by the theory are not istateed with, the 
steadiness of the layer will be assured. 


The author acknowledges gratefully the editorial assistance in the 
preparation of this paper furnished by George W. Patterson, and 
also the valuable suggestions offered by Galen B. Schubauer. 
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1 to 5.6 at intervals of 0.001. 
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x= {0(.001)7(.01)50(.1)300(1)2,000(10)10,000;12D} 
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5 pages with cover. 25 cents, 
MT19. On the Function H (m, a, x)=exp(—ix) F (m+-1—ia, 2m-+-2; ix); with table of the 
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(Reprinted from J. Math. Phys., December 1942.). 20 pages, with cover, 25 cents. 


MT20. Table of integrals f Jo(t)dt and ligvors 


Values of the two integrals are given for x==0(.01)10 to 10 decimal places. (Reprinted from 
J. Math. and Phys., May 1943.) 12 pages, with cover, 25 cents. 


MT21. Table of Jio(x)= f, “1, and Related Functions: 
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differences of F(x). 


with even central 

Table II: Jig(x) to 10 decimal places, for x==3(.1)10(1)22 with even central differences up to 
x=100. 

Table III: “Reduced” derivatives of F(x) for x= 10(1)21 and n==0(1)13, to 12 decimal places. 

(Reprinted from J. Math. Phys., June 1943.) 7 pages, with cover, 25 cents. 

MT22. Table of Coefficients in Numerical Integration Formulae: 

The values of B‘”) ,(1)/n! and B‘™,, /n! where B‘™ ,,(1) denotes the n*® Bernoulli polynomial of 
the nth order for x=1 and B”,, denotes the n*® Bernoulli number:of the nt order, were com- 
puted for n=1, 2,... 20. The quantities B‘” ,(1)/n! are required in the Laplace formula of | 
numerical integration “employing forward differences, as well as in the Gregory formula. The 
quantities B‘”) ,,/n! are used in the Laplace formula employing backward differences. 

(Reprinted from J. Math. Phys., June 1943.) 2 pages, with cover, 25 cents. 

MT23. Taste or Fourter CozrricienTs. 

Whenever ¢ (x) is a known polynomial whose degree does not exceed 10, the present table 

of the functions 


1 1 
shn)= f x* sin nex dx and Chn= f x* cos nx dx to 10D(1==k=10, 1 =—n=100), will 
facilitate the evaluation of the first hundred Fourier Coefficients. 
Reprinted from J. Math. Phys. Sept., 143.) 11 pages, with cover, 25 cents. 


MT24. Coerricirents ror Numericat DirrerenTIATION Wira Centrat DirreReNnces. 
Coefficients are given for derivatives as far.as the 52d. For the first 30 derivatives, exact 
values are given for coefficients of the first 30 differences, and also exact values are given for some 
coeflicients of differences beyond the 30th. For the other coefficients, values are given to 18 


t figures 
(Reprinted fro from J. Math. Phys., Sept. 1943) 21 pages, with cover, 25 cents. 


Payment is required in advance. Make remittance payable to the “National 
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